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Abstract 

We  present  theoretical  results  on  the  stability  properties  of  a  Blasius  boundary  layer  perturbed  by 
Klebanoff  distortions  with  a  relatively  long  spanwise  scale.  Even  relatively  weak  Klebanoff  modes  can 
alter  the  near- wall  curvature  of  the  underlying  flow  by  0(1)  and,  hence,  introduce  linear  instabilities  with 
larger  characteristic  growth  rates  and  frequencies  than  those  of  the  Tollmien-Schlichting  waves  in  an 
unpertubed  Blasius  flow.  A  localised  distortion  supports  both  sinuous  and  varicose  modes  of  instability, 
with  the  growth  rates  of  the  sinuous  modes  being  likely  to  be  larger,  in  general.  Overall,  the  instability  is 
intermittent  in  time  and  localised  in  space,  being  confined  to  a  finite  part  of  the  Klebanoff  mode  cycle  and 
to  a  specific  window(s)  along  the  streamwise  direction.  A  spanwise  periodic  distortion  supports  spatially 
quasi-periodic  modes  (via  the  parametric  resonance  mechanism),  which  may  be  viewed  as  modified  T-S 
waves  with  excess  growth  rates  when  the  Klebanoff  modes  are  weak.  In  spite  of  the  simplifications 
involved  in  this  theory,  its  predictions  appear  qualitatively  consistent  with  some  of  the  unusual 
characteristics  of  the  high-frequency  wavepackets  observed  during  previous  experiments.  The  nonlinear 
development  of  a  localised  sinuous  mode  is  followed  across  a  sequence  of  asymptotic  regimes  using  the 
non-equlibrium  critical-layer  theory. 
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LINEAR  AND  NONLINEAR  INSTABILITIES  OF  BLASIUS  BOUNDARY 
LAYER  PERTURBED  BY  STREAMWISE  VORTICES 
PART  II:  INTERMITTENT  INSTABILITY  INDUCED  BY 
LONG- WAVELENGTH  KLEBANOFF  MODES 

XUESONG  WU*  AND  MEELAN  CHOUDHARlt 

Abstract.  This  paper  presents  theoretical  results  on  the  stability  properties  of  a  Blasius 
boundary  layer  perturbed  by  Klebanoff  modes.  The  latter  are  taken  to  be  the  signature  of  a 
low-frequency  three-dimensional  convected  gust,  which  may  be  either  isolated  or  periodic  along 
the  spanwise  direction.  They  can  be  treated  as  a  distortion  to  the  basic  state  because  of  their 
extremely  low  frequency  and  long  streamwise  wavelength.  It  is  found  that  even  relatively  weak 
Klebanoff  modes  may  alter  the  curvature  of  the  underlying  mean  flow  by  0(1)  in  a  near- wall  region. 
The  resulting  perturbed  flow  may  support  linear  instability  modes  with  a  streamwise  wavelength 
much  shorter  than  the  spanwise  length  scale  of  the  distortion. 

When  the  amplitude  of  the  Klebanoff  modes  exceeds  a  threshold  range,  the  instability  modes 
become  predominantly  inviscid,  and  their  growth  rates  and  characteristic  frequencies  are  much 
higher  than  those  of  the  Tollmien-Schlichting  (T-S)  waves  in  an  unperturbed  Blasius  flow.  A 
localised  distortion  supports  both  sinuous  and  varicose  modes  of  instability,  with  the  sinuous  modes 
being  generally  more  unstable  than  the  varicose  modes.  Overall,  the  instability  is  intermittent  in 
time  and  localised  in  space,  occurring  only  in  certain  phases  of  the  modulation  cycle  and  within 
a  specific  window(s)  along  the  streamwise  direction.  In  particular,  the  dominant  sinuous  modes 
appear  only  during  the  phase  in  which  a  low-speed  streak  exists.  A  periodic  distortion  supports 
spatially  quasi-periodic  modes  through  a  parametric  resonance  mechanism. 

The  nonlinear  development  of  a  localised  sinuous  mode  is  followed  using  nonequilibrium  critical- 
layer  theory.  During  the  first  stage  of  nonlinear  evolution,  the  amplitude  function  is  governed  by  an 
integro-ordinary-differential  equation,  the  solution  of  which  develops  a  finite-distance  singularity  if 
the  viscous  effects  are  not  too  large.  In  the  vicinity  of  this  singularity,  the  disturbance  amplifies  over 
a  much  shorter  streamwise  scale.  Its  amplitude  is  now  governed  by  an  integro-partial-differential 
equation  instead,  which  describes  the  (simultaneous)  spatial  and  nonlinear  modulation  of  the  in¬ 
stability  mode.  The  solution  of  this  equation  appears  to  develop  yet  another  singularity,  indicating 
the  onset  of  fully  nonlinear  yet  inviscid  dynamics.  We  suggest  that  the  temporally  intermittent  and 
spatially  localised  instability  and  its  subsequent  nonlinear  development  lead  to  patches  of  streak 
oscillations.  The  apparent  convection  velocity  of  such  patches  is  estimated  to  be  about  0.82  times 
the  free-stream  velocity. 
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This  Research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract  No. 
NAS  1-97046  while  in  residence  at  ICASE,  NASA  Langley  Research  Center,  Hampton,  VA  23681-2199. 

f  Computational  Modeling  and  Simulation  Branch,  NASA  Langley  Research  Center,  Hampton,  VA 
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1.  Introduction.  It  is  known  that  low-frequency  components  of  three-dimensional  vortical 
disturbances  in  the  free  stream  can  be  entrained  into  the  boundary  layer  due  to  the  nonparallel 
flow  effect,  producing  significant  distortion  in  form  of  alternate  thickening  and  thinning  of  the 
layer  along  the  spanwise  direction.  This  observation  goes  back  to  Dryden  [9]  and  Taylor  [31]  who, 
in  fact,  suggested  that  the  entrained  vortex  motion,  rather  than  the  Tollmien-Schlichting  (T-S) 
instability,  was  the  leading  cause  behind  transition  to  turbulence.  The  dispute  continued  until  the 
experiments  of  Schubauer  &  Skramstad  [29],  which  fully  validated  the  work  of  Tollmien  [33]  and 
Schlichting  [28]  by  suppressing  the  effects  of  free-stream  turbulence.  Since  then,  most  studies  have 
focused  on  transition  under  low-amplitude  disturbance  environments,  for  which  the  T-S  modes 
provide  the  starting  point  for  understanding  the  overall  transition  process. 

Motivated  by  its  relevance  to  turbomachinery,  however,  there  has  also  been  a  considerable 
amount  of  work  on  transition  at  moderate  to  high  levels  of  free-stream  turbulence.  This  has 
renewed  the  interest  in  the  findings  of  Dryden  [9]  and  Taylor  [31].  Recent  experimental  studies 
(see  e.g.  Kendall  [18],  Westin  et  al.  [37],  Matsubara  &  Alfredsson  [26],  and  references  therein)  show 
that  the  boundary  layer  filters  out  the  high-frequency  components  of  free-stream  turbulence,  while 
amplifying  the  low-frequency  parts  of  the  signature.  The  distortion  within  the  boundary  layer  is 
dominated  by  the  streamwise  velocity  fluctuations,  which  are  manifested  in  the  form  of  longitudinal 
vortices  or  streaks.  In  recognition  of  the  contribution  by  Klebanoff  [22],  Kendall  [18]  referred  to 
these  low-frequency  streaks  as  the  Klebanoff  modes,  and  that  name  has  been  widely  adopted  in 
the  literature  since  then.  However,  because  the  Klebanoff  modes  do  not  actually  represent  the 
eigenmodes  of  the  flow,  we  shall  refer  to  them  as  Klebanoff  distortions  or  fluctuations  in  order  to 
avoid  potential  confusion  when  genuine  instability  modes  are  being  discussed. 

Since  Klebanoff  fluctuations  are  of  low  frequency,  one  might  be  inclined  to  believe  that  an 
appropriate  steady  perturbation  may  capture  the  essential  physics  involved.  Indeed,  a  steady 
perturbation  in  the  oncoming  flow  can  induce  a  boundary-layer  distortion  that  is  similar  to  the 
Klebanoff  motions  (Bradshaw  [5]).  Initiated  by  Bradshaw’s  observations,  Crow  [8]  calculated  the 
resulting  distortion  using  linearized  boundary-layer  equations  and  was  able  to  describe  the  alter¬ 
nate  thickening  and  thinning  of  the  boundary  layer  near  the  leading  edge  as  noted  in  the  earlier 
experiments  by  Taylor  [31]  and  Bradshaw  [5].  The  boundary-layer  response  to  small-amplitude 
unsteady  vortical  disturbances  was  calculated  by  Gulyaev  et  al.  [16]  and  Choudhari  [6]  using  the 
unsteady  form  of  linearized  boundary-layer  equations.  They  showed  that  the  low-frequency  distur¬ 
bances  induce  a  large  streamwise  velocity  fluctuation  in  the  boundary-layer  region,  which  exhibits 
a  significant  growth  in  the  downstream  direction. 

The  boundary-layer  approximation  is  valid  sufficiently  close  to  the  leading  edge,  where  the 
thickness  of  the  boundary  layer  is  smaller  than  the  spanwise  length  scale  of  the  perturbation.  Due 
to  a  continued  growth  of  the  boundary-layer  thickness  (as  well  as  of  the  perturbation  amplitude), 
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however,  cross-flow  ellipticity  (and  disturbance  nonlinearity)  must  become  significant  sufficiently 
downstream  of  this  region.  The  subsequent  evolution  of  the  boundary-layer  signature  of  free-stream 
disturbances  is  described  by  the  boundary-region  equations  (Goldstein  &  Leib  [13];  Wundrow  & 
Goldstein  [43];  and  Leib,  Wundrow  &  Goldstein  [24]).  Wundrow  &  Goldstein  [43]  considered  the 
generic  case  for  purely  steady  perturbations,  such  that  the  nonlinearity  and  cross-flow  ellipticity 
come  into  play  simultaneously.  In  this  case,  the  distorted  flow  differs  from  the  original  Blasius 
profile  by  0(1)  and,  moreover,  becomes  inflectional  in  the  main  part  of  the  boundary  layer.  It  can 
therefore  support  a  new  and  stronger  class  of  inviscid  Rayleigh  instabilities,  which  may  lead  to 
the  onset  of  turbulence  via  an  alternative  route  to  the  T-S  modes,  namely,  the  by-pass  transition. 
Leib  et  al  [24]  used  the  linearized  form  of  the  unsteady  boundary-region  equations  to  describe  the 
low-frequency  Klebanoff  fluctuation  further  downstream  of  the  leading  edge.  By  representing  the 
broad-band  free-stream  turbulence  as  a  superposition  of  Fourier  modes,  Leib  et  al.  [24]  calculated 
the  streamwise  evolution  of  the  root-mean-square  of  the  fluctuations  within  the  boundary  layer, 
which  compared  favourably  with  the  relevant  experimental  data. 

As  indicated  above,  experiments  have  provided  fairly  complete  quantitative  data  about  the 
characteristics  of  Klebanoff  fluctuations  themselves,  and  considerable  theoretical  progress  has  also 
been  made  in  terms  of  characterising/predicting  them.  However,  the  transition  process  in  the 
presence  of  these  fluctuations  remains  poorly  understood.  In  particular,  it  is  not  clear  if  the  T-S 
waves  continue  to  play  an  important  part  in  transition  or  if  the  latter  is  initiated  by  some  other 
instabilities,  e.g.  the  instability  of  the  streaks. 

Direct  laboratory  investigations  of  the  transition  process  in  the  presence  of  Klebanoff  fluctua¬ 
tions  have  been  made  by  a  number  of  investigators.  At  moderate  levels  of  free-stream  turbulence, 
Arnal  &  Juillen  [2]  and  Kendall  [19]  have  observed  an  intermittent  appearance  of  wavepackets 
inside  the  boundary  layer.  While  the  exact  origin  of  the  wavepackets  was  unclear,  they  appar¬ 
ently  resembled  T-S  waves.  A  series  of  landmark  experiments  conducted  by  Kendall  has  revealed 
some  unusual  characteristics  of  these  wavepackets.  First,  they  appear  only  when  the  free-stream 
turbulence  level  exceeds  a  threshold  of  about  0.1%.  Second,  the  growth  rates  of  these  packets 
are  considerably  larger  than  those  of  the  T-S  waves,  being  also  dependent  on  the  turbulence  level 
(Kendall  [20],  [21]).  Third,  the  frequencies  of  these  waves  are  also  higher  than  those  of  T-S  instabil¬ 
ity.  Finally,  these  packets  are  relatively  confined  in  the  spanwise  direction,  spreading  rather  slowly 
as  they  propagate  downstream.  There  has  been  no  physical  explanation  for  the  above  features  thus 
far,  except  for  the  conjecture  by  Goldstein  &  Wundrow  [15],  who  attributed  the  intermittency  of 
the  wavepackets  and  the  nonlinear  dependence  of  their  amplitude  on  the  turbulence  level  to  ran¬ 
dom  receptivity  via  nonparallelism  of  the  mean  flow  (due  to  either  a  rapid  boundary-layer  growth 
near  the  leading  edge  and/or  geometric  perturbations  in  the  region  downstream).  An  alternative 
hypothesis  (which  does  not  preclude  the  scenario  proposed  by  Goldstein  &  Wundrow)  is  suggested 
in  the  present  study,  which  attempts  to  shed  further  light  on  some  of  the  unusual  characteristics 
of  the  high-frequency  wavepackets  noted  above. 

Rather  than  studying  naturally  occurring  wavepackets,  Watmuff  [34]  used  a  harmonic  point 
source  to  generate  spanwise  localised  wavetrains  in  a  controlled  fashion.  Even  a  weak  Klebanoff 
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perturbation  was  found  to  severely  distort  the  wavetrain,  such  that  a  comparison  with  calculations 
that  did  not  account  for  the  Klebanoff  distortions  was  almost  meaningless.  On  the  other  hand, 
however,  the  artificially  excited  wavetrain  did  not  exhibit  the  excess  growth  as  observed  by  Kendall 
for  the  naturally  occurring  wavepackets.  Watmuff  attributed  this  difference  to  the  fact  that  the 
Klebanoff  fluctuations  were  relatively  weak  in  his  experiments. 

Boiko  et  al.  [4]  studied  the  development  of  an  artificially  excited  planar  T-S  wave  in  a  boundary 
layer  subjected  to  1.5%  free-stream  turbulence.  As  expected,  the  planar  wave  front  was  deformed 
by  the  Klebanoff  fluctuations;  however,  the  growth  rate  of  the  T-S  wave  was  reduced  (relative 
to  that  in  an  unperturbed  Blasius  flow),  which  suggests  that  the  Klebanoff  fluctuation  now  had 
a  stabilising  effect  on  the  unstable  modes  at  higher  frequencies.  The  ultimate  breakdown  of  the 
laminar  boundary  layer  was,  however,  still  associated  with  the  nonlinear  development  of  the  T-S 
wave.  These  observations  present  a  rather  paradoxical  situation,  since  it  is  well  known  that  higher 
intensity  free-stream  turbulence  generally  hastens  the  laminar-turbulent  transition.  A  similar  prob¬ 
lem  was  recently  investigated  by  Watmuff  [36],  albeit  in  a  different  experimental  setting.  He  also 
observed  the  distortion  of  the  wavefront,  but  found  that  the  growth  rate  was  hardly  affected  by 
the  Klebanoff  distortions.  Nevertheless,  transition  in  his  experiments  did  occur  slightly  earlier 
than  when  the  Klebanoff  fluctuation  was  absent.  Finally,  Bakchinov  et  al.  [3]  have  investigated 
the  interaction  between  a  planar  T-S  wave  of  relatively  high  frequency  and  a  localised  disturbance 
introduced  in  the  free-stream  region.  When  generated  separately,  each  of  these  disturbances  was 
found  to  decay.  But  when  introduced  simultaneously,  Bakchinov  et  al.  observed  enhanced  ampli¬ 
fication  of  a  broad-band  low-frequency  oblique  waves,  presumably  due  to  a  nonlinear  interaction 
between  the  two  artificially  introduced  disturbances.  The  nature  of  this  interaction  has  not  been 
explained  as  yet. 

In  addition  to  investigating  the  role  of  conventional  T-S  waves  during  transition  under  moderate 
levels  of  free-stream  turbulence,  it  is  also  important  to  examine  alternative  instability  mechanisms. 
Streak  instability,  in  particular,  has  attracted  a  lot  of  attention  in  recent  years.  Matsubara, 
Bakchinov  &  Alfredsson  [25]  reported  that  streaks  or  Klebanoff  modes  are  unstable,  which  can 
lead  to  a  meandering  and  oscillation  of  the  streaks  and  an  eventual  breakdown  into  turbulent 
spots.  Unfortunately,  the  random  nature  of  the  streaks  implies  that  a  detailed  quantitative  study 
of  their  stability  would  be  rather  difficult.  Therefore,  some  researchers  have  chosen  to  generate 
steady  spanwise-dependent  flow  in  a  controlled  manner,  typically  through  surface  suction/blowing 
or  surface  roughness.  These  artificially  created  distortions  are  akin  to  Klebanoff  modes  in  some 
respects  and,  hence,  these  kinds  of  studies  may  shed  useful  light  on  the  streak  instability.  For  a 
survey  of  the  main  findings  from  these  studies,  the  reader  is  referred  to  the  Part  I  of  this  paper 
(Wu  &  Luo  [40]). 

In  an  effort  to  understand  streak  breakdown,  Andersson  et  al.  [1]  modelled  the  streak  structure 
by  a  steady,  spanwise-periodic  distortion  to  a  Blasius  boundary  layer.  An  inviscid  stability  analysis 
based  on  Floquet  theory  suggested  that  the  streaks  become  unstable  only  when  the  amplitude  of  the 
associated  streamwise  velocity  perturbation  exceeds  approximately  26%  of  the  free-stream  velocity. 
This  estimate  is  perhaps  too  high  to  be  representative  of  typical  Klebanoff  distortion  in  natural 
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disturbance  environments.  Based  on  the  present  work,  it  appears  that  a  partial  explanation  for 
this  discrepancy  may  well  be  related  to  the  unsteadiness  of  the  Klebanoff  distortion,  which  (in  spite 
of  the  low  frequencies  of  the  distortion)  is  found  to  exert  a  significant  effect  on  the  high-frequency 
secondary  instabilities  (see  §2.2). 

A  direct  numerical  simulation  of  the  transition  initiated  by  high  levels  of  free-stream  turbulence 
was  first  performed  by  Rai  &  Moin  [27].  Their  simulations  involved  three  sequential  streamwise 
domains  in  order  to  cope  with  the  different  resolution  requirements  associated  with  the  three 
dominant  physical  processes  relevant  to  this  problem,  viz.  (a)  interaction  of  free-stream  turbulence 
with  the  leading  edge,  (b)  subsequent  development  of  disturbance  within  the  boundary  layer  plus 
the  early  transition  stage,  and  (c)  the  late  stages  of  transition  plus  the  fully  turbulent  state.  This 
treatment  allowed  for  a  realistic  description  of  the  practical  situation  while  avoiding  the  excessive 
grid  count  that  would  otherwise  have  been  required. 

Jacobs  &  Durbin  [23]  carried  out  a  direct  numerical  simulation  of  the  bypass  transition  caused 
by  strong  free-stream  disturbances,  which  were  represented  by  a  superposition  of  modes  from  the 
continuous  spectrum  of  the  Orr-Sommerfeld  equation.  Unlike  Rai  &  Moin  [27],  the  perturbation 
was  imposed  at  some  large  distance  downstream  from  the  leading  edge.  Due  to  the  reduced  size 
of  the  simulation  domain  and  the  increased  computational  power  over  the  years,  Jacobs  &  Durbin 
were  able  to  perform  the  simulation  on  a  finer  mesh  than  what  was  available  to  Rai  &  Moin 
[27].  It  was  found  that  the  boundary-layer  response  was  indeed  dominated  by  streamwise  streaks. 
However,  these  streaks  appeared  to  be  fairly  stable  and  it  was  only  after  they  had  lifted  up  to  the 
outer  part  of  the  boundary  layer  to  form  a  ‘backward  jet’  that  the  breakdown  to  turbulent  spots 
occurred.  Their  simulation  also  revealed  that  the  breakdown  was  local,  being  determined  by  a 
combination  of  the  strength  and  the  length  scale  of  the  local  distortion;  the  collective  instability 
of  the  Floquet  type  was  not  observed. 

In  this  paper,  we  investigate  the  instability  of  a  Blasius  boundary  layer  perturbed  by  Klebanoff 
distortions.  Our  main  interest  will  be  in  small-amplitude  distortions,  partly  because  they  allow  the 
phenomenon  to  be  studied  on  an  analytical  basis  but  also  because  the  amplitude  levels  measured 
in  most  experiments  do  not  exceed  10-15%  of  the  free-stream  velocity.  One  significant  assumption 
underlying  the  present  work  involves  the  spanwise  length  scale  of  the  Klebanoff  distortion,  which  is 
assumed  to  be  larger  than  the  boundary-layer  thickness.  As  in  Part  I,  we  address  two  main  issues: 
(a)  how  the  Tollmien-Schlichting  instability,  which  operates  in  the  absence  of  any  distortion,  is 
modified  by  a  weak  Klebanoff  fluctuation,  and  (b)  whether  a  weak  distortion  can  induce  an  inviscid 
instability  which  is  absent  from  the  unperturbed  boundary  layer.  An  asymptotic  approach  based 
on  the  high-Reynolds-number  assumption  is  employed  to  describe  both  the  Klebanoff  fluctuation 
and  the  instability  of  the  perturbed  flow  in  a  systematic  and  consistent  manner. 

In  §2,  we  formulate  the  problem  in  an  asymptotic  framework.  A  three-dimensional  vortical 
disturbance  is  prescribed  in  the  oncoming  flow.  The  spanwise  distribution  of  the  perturbation 
is  allowed  to  be  either  localised  or  periodic.  The  relevant  inviscid  solution,  valid  outside  of  the 
boundary  layer,  is  obtained  first.  For  the  scalings  adopted  in  this  paper,  the  boundary-layer 
signature  of  this  inviscid  perturbation,  i.e.  the  Klebanoff  distortion,  is  governed  by  the  linearized 


5 


unsteady  boundary-layer  equations.  In  §2.2,  we  present  the  details  of  the  scaling  arguments  which 
determine  both  the  characteristic  wavelength  of  the  instability  modes  as  well  as  the  required 
magnitude  of  the  Klebanoff  distortion.  Depending  on  the  relative  values  of  the  streamwise  and 
spanwise  length  scale  of  the  instability,  two  distinct  types  of  instability  modes  may  be  considered. 
The  first  class  of  modes  corresponds  to  the  fully  three-dimensional  (i.e.  oblique)  modes  that  were 
identified  by  Goldstein  and  Wundrow  [14]  in  the  context  of  purely  stationary  distortions  to  a 
Blasius  flow.  The  other  set  of  modes  is  primarily  two-dimensional,  i.e.  it  has  significantly  shorter 
wavelengths  in  the  streamwise  direction  than  the  spanwise  length  scale  involved.  The  latter  scale 
is,  of  course,  the  same  as  that  of  the  distortion  for  both  classes  of  modes.  The  Goldstein- Wundrow 
(GW)  modes  are  based  on  a  parametric  resonance  mechanism,  and  therefore  exist  only  for  global 
distortions  that  are  periodic  in  the  spanwise  direction. 

The  linear  instability  of  the  modified  flow  is  analysed  in  §3.  Due  to  our  interest  in  high- 
frequency  wavepackets  that  are  localised  along  the  spanwise  direction,  the  case  of  localised  distor¬ 
tions  (and,  therefore,  the  second  class  of  instabilities  mentioned  above)  is  emphasised  in  this  paper. 
We  show  that  the  spanwise  mode  shape  of  the  instability  is  governed  by  a  Schrodinger  equation, 
which  describes  the  interaction  between  the  instability  mode  and  the  Klebanoff  distortion  within 
the  bulk  of  the  boundary  layer.  The  potential  function  in  this  Schrodinger  equation  is  proportional 
to  the  local  wall-shear  perturbation  induced  by  the  the  Klebanoff  distortion.  The  growth  rate  of 
the  instability  mode  is  determined  at  a  higher  order  in  the  asymptotic  analysis,  by  the  curvature  of 
the  distorted  flow  at  the  critical  level  of  the  instability  mode.  For  a  localised  distortion,  solutions 
representing  sinuous  and  varicose  modes  are  presented,  and  their  relative  importance  is  assessed. 
Further  calculations  are  carried  out  for  a  spanwise  localised  Klebanoff  fluctuation  to  demonstrate 
the  intermittent  and  local  nature  of  the  instability. 

The  nonlinear  development  of  a  sinuous  instability  mode  through  a  sequence  of  distinct  asymp¬ 
totic  sub-regimes  is  considered  in  §4.  In  §5,  we  discuss  how  the  preceding  analysis  can  be  applied 
to  the  case  of  spanwise-periodic  Klebanoff  distortions.  The  main  findings  from  this  work  and  their 
physical  implications  are  summarized  in  §6. 

2.  Formulation.  Consider  the  two-dimensional  incompressible  boundary  layer  due  to  a  uni¬ 
form  flow  with  velocity  UD 0  past  a  semi-infinite  flat  plate.  Superimposed  on  the  incoming  stream 
is  a  small-amplitude,  three-dimensional  vortical  disturbance  (i.e.  a  ‘gust’)  that  is  assumed  to  be 
advected  at  the  free-stream  speed.  For  simplicity,  we  also  assume  that  the  free-stream  disturbance 
is  harmonic  in  time,  with  a  frequency  of  kiUoo/A  where  k\  denotes  the  nondimensional  frequency 
parameter  and  A  represents  the  dimensional  length  scale  of  the  gust  in  the  spanwise  direction. 

The  flow  is  described  in  the  Cartesian  coordinate  system  (x,y,  Z),  which  has  its  origin  at  the 
plate  leading  edge.  Here,  x,  y  and  Z  denote  the  streamwise,  normal,  and  spanwise  coordinates 
nondimensionalized  with  respect  to  A.  The  time  variable  t  is  normalized  by  A /JJoq  and  the  velocity 
components  (u,  v,  w)  and  pressure  p  are  normalized  by  Uoq  and  pU^  respectively,  where  p  denotes 
the  fluid  density.  The  Reynolds  number 
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is  assumed  to  be  a  large  parameter  throughout  this  analysis  (i.e.  R\  1). 

The  streamwise  and  normal  velocity  components  of  the  Blasius  flow  are  given  by 

(UB,VB)  =  {f'(ti),  (2 xRA)-i(r,F' -F)}  , 

where  F  satisfies  the  Blasius  equation 

(2.1)  F'"  +  FF"  =  0  ,  F(0)  =  F'( 0)  =  0  ,  F'( oo)  1  ,  with  r\  =  R\^=  . 

V2x 

2.1.  Flow  distortion  induced  by  KlebanofF  modes.  The  boundary-layer  response  to  a 
three-dimensional  convected  gust  was  analyzed  by  Gulyaev  et  al  [16],  Choudhari[6]  and  Leib  et 
al  [24] .  Similar  to  their  work,  the  velocity  field  of  the  disturbance  superimposed  on  the  oncoming 
flow  has  the  form 

(2.2)  Uoo  =  eD(u00B'(Z),v00B'(Z),w00B(Z ))  ? 

where  eo  represents  the  gust  amplitude  (cd  <  1),  and  k\  and  fe  denote  the  streamwise  and 
transverse  wavenumbers,  respectively.  As  shown  in  the  above  studies,  the  amplification  of  the 
gust  signature  within  the  boundary  layer  is  directly  proportional  to  the  ratio  of  the  spanwise 
and  streamwise  wavenumbers,  respectively.  Accordingly,  we  assume  that  Aq,  <  1,  which  also 
amounts  to  Uoq,  Voq  (=  0(1)).  For  the  most  part,  our  interest  will  be  confined  to  the  range 

of  streamwise  locations  that  corresponds  to 

(2.3)  £  <  Ra  , 

so  that  the  spanwise  length  scale  of  the  disturbance  is  significantly  greater  than  the  local  thickness 
of  the  boundary-layer,  and 

(2.4)  x  <C  Cp1 ,  , 

which  allows  one  to  suppress  the  effects  of  disturbance  nonlinearity  as  the  gust  is  advected  with 
the  free  stream  (Leib  et  al.  [24]).  We  observe  that  the  measured  properties  of  Klebanoff  modes 
are  closely  tied  to  the  wavenumber-frequency  spectrum  of  turbulence  in  the  upstream  flow  and, 
therefore,  will  vary  from  facility  to  facility  as  well  as  between  wind  tunnels  and  the  flight  environ¬ 
ment.  Therefore,  the  accuracy  of  the  “long  spanwise- wavelength”  assumption  will  accordingly  vary 
from  one  disturbance  environment  to  another.  For  example,  as  pointed  out  by  Leib  et  al  [24],  the 
measured  length  scale  of  Klebanoff  fluctuations  was  nearly  five  times  larger  than  the  boundary- 
layer  thickness  in  the  original  experiments  of  Klebanoff  [22];  however,  the  scale  was  closer  to  the 
boundary-layer  thickness  in  the  experiments  by  Kendall  [18]  and  Westin  et  al  [37].  In  the  present 
context,  of  course,  the  assumptions  (2.3)  and  (2.4)  have  been  motivated  primarily  by  the  resulting 
simplicity  of  the  analysis. 

Note  that  we  have  allowed  for  a  general  spanwise  dependence  of  the  gust  via  the  arbitrary 
function  B(Z).  This  generalization  was  prompted  by  the  fact  that,  in  practice,  the  free-stream 
disturbances  are  stochastic  in  nature,  having  a  finite  correlation  length  along  the  spanwise  direction. 
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Any  given  realization  may,  therefore,  resemble  a  finite-extent  perturbation  rather  than  a  spanwise 
periodic  one.  A  similarly  general  dependence  of  the  free-stream  disturbance,  either  on  x  (or, 
equivalently,  on  t  due  to  the  convective  phase)  and/or  on  y ,  may  also  be  allowed  for  in  this 
framework  via  a  Fourier  superposition.  However,  for  simplicity,  we  will  proceed  on  the  basis  of 
the  harmonic  assumption  (2.2)  and,  later,  indicate  how  the  final  results  may  be  generalized  to  an 
arbitrary  form  of  convected  free-stream  vorticity. 

The  inviscid  solution  (valid  outside  of  the  main  boundary  layer)  for  k\X  —  0(1)  can  be  written 
as  (Leib  et  al.  [24]) 

(2.5)  up  =  +  y  0] 

where  denotes  the  velocity  potential  and,  to  the  required  order,  it  satisfies  the  Laplace  equation 

(2.6)  V20  =  0  ■ 


The  boundary  conditions  are  given  by 


(2.7) 


V0  — y  0  as  y  oc  ; 

<j>(x.  0)  =  0  (x  <  0),  and  0y(x,O )  =  —v00B,(Z)eikl^x~t^  (x  >  0) 


The  full  solution  to  the  boundary-value  problem  (2.6)-(2.7)  can  be  found  by  the  standard  Wiener- 
Hopf  technique.  However,  for  the  purpose  of  stability  analysis,  we  confine  ourselves  to  the  region 
x  >  1,  such  that  k\x  —  0(1)  while  (2.3)  is  still  satisfied.  The  inviscid  solution  under  these 
conditions  can  be  obtained  by  neglecting  the  x-derivative  term  in  (2.6)  and  solving  the  resultant 
two-dimensional  Laplace  equation  in  the  half  space  y  >  0.  This  yields  the  following  solution  for 
the  slip  velocity  components  in  streamwise  and  spanwise  directions: 


(2.8) 


Ua 


Wb 


d_r 

dzj_c 


m 

C-z 


dC- 


The  signature  of  the  above  inviscid  solution  within  the  boundary  layer  corresponds  to  a  small 
perturbation  to  the  Blasius  flow  and,  to  the  leading  order  of  approximation,  is  given  by 

(2.9)  uD=eD[U,(^)1W,W]e~i*+c.c.  +  ...  , 

Ua 

where  we  have  put 


x  =  Aqx,  t  =  kit  . 

Irrespective  of  the  value  of  the  transverse  wavenumber  fe  or  the  spanwise  distribution  B(Z ),  the 
solution  for  the  boundary-layer  distortion  can  be  split  into  a  sum  of  two  parts: 

(2.10)  (U,V,W)  =  [-^Zlu,-^Zlv,ws(Z)w]  +us[Ui,V!,0]  +  ...  . 

Since  k\  <  1  for  the  low-frequency  Klebanoff  fluctuation,  the  second  part  is  asymptotically  smaller 
than  the  first  and,  therefore,  does  not  influence  the  stability  of  the  distorted  flow  to  the  leading 
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1.5 


Fig.  2.1.  Variation  of  U"'(0,x)  (solid  line:  real  part,  dashed  line:  imaginary  part)  and  \Uf(0,x)\  (—•  —  •—) 
with  x. 


order.  Accordingly,  only  the  first  part  of  the  Klebanoff  distortion  in  (2.10)  needs  to  be  considered 
in  the  present  work.  Equation  (2.10)  also  shows  that  the  streamwise  velocity  perturbation  within 
the  boundary  layer  is  much  larger  in  amplitude  than  the  free-stream  disturbance.  The  functions 
(U,V,W)  are  governed  by  the  linearized  unsteady  boundary-layer  equations  (Leib  et  al.  [24]) 


(2.11) 

(2.12) 

(2.13) 


-i  U  +  F'^r 


dU_rLdU  W 

dx  2x  dr}  dr)  ’ 

_Z^..hHu+f,,v  =  -— 

2x  drj  2x  2x  dr]2 

-  ,dW  F  dW  1  d2W 

-\W  - h  F - = - . 

dx  2x  dr}  2x  dr}2 


where  r}  is  defined  by  (2.1).  The  appropriate  boundary  conditions  are 


(2.14) 


U  =  V  =  W  =  0  at  r}  =  0  , 

£7  — ^  0,  as  r}  oc. 


In  the  upstream  limit  (x  0),  the  flow  becomes  quasi-steady  and  its  solution  matches  with  that 
obtained  by  Crow  [8], 

(2.15)  U  ^  \xrtF"  ,  V  ^hr)2F"  -3r)F'  -F)  ,  W  ->■  F'  . 

Z  ti 

On  the  other  hand,  as  x  oc,  the  perturbations  move  away  from  the  wall  and  towards  the  edge 
of  the  mean  boundary  layer  (Choudhari  [6],  Leib  et  al.  [24]).  Because  the  maximum  distortion  to 
the  boundary-layer  flow  is  confined  to  the  region  x  =  0(1),  this  is  the  region  of  interest  from  the 
standpoint  of  investigating  potentially  radical  changes  to  the  linear  instability  of  the  flow. 

We  will  show  in  §3  that  the  complete  profiles  of  the  Klebanoff  distortion  are  not  required  to 
determine  the  leading  order  instability  characteristics  of  the  perturbed  flow;  instead,  the  charac¬ 
teristics  are  completely  determined  by  local  values  of  the  slope  U'( 0,x)  and  the  torsion  U,ff( 0,x) 
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of  the  streamwise  velocity  fluctuation  near  the  surface.  In  fact,  the  wall  torsion  itself  is  related  to 
the  wall-shear  perturbation  by  the  simple  relation 

(2.16)  U'"(0,x)  =  -2ixU'(0,x)  , 

which  follows  from  a  differentiation  of  (2.11)  with  respect  to  rj  and  setting  rj  =  0.  The  importance 
of  the  near- wall  features  of  the  gust  signature  is  easily  understood  by  recalling  that  the  viscous 
instability  modes  of  the  unperturbed  Blasius  flow  are  controlled  by  the  mean  value  of  the  skin 
friction  parameter  and,  as  noted  by  Goldstein  &  Wundrow  [14],  any  inflectional  behaviour  of  the 
perturbed  base  flow  is  also  confined  to  a  narrow  vicinity  of  the  surface  (since  the  amplitude  of 
the  Klebanoff  distortion  is  assumed  to  be  small  in  comparison  with  the  free-stream  speed).  The 
streamwise  variation  in  the  amplitudes  of  the  wall-shear  and  wall-torsion  fluctuations  due  to  the 
gust  is  displayed  in  figure  2.1,  which  underscores  the  previously  made  observation  that  the  region 
x  =  0(1)  accounts  for  the  range  of  locations  (or,  equivalently,  the  range  of  frequencies  at  a  given 
streamwise  location)  wherein  the  boundary-layer  flow  is  distorted  the  most. 

2.2.  Scaling  arguments.  As  discussed  in  §2.1,  the  Klebanoff  distortion  in  the  boundary 
layer  is  concentrated  in  the  streamwise  region  where  x  =  0(1),  i.e.  at  a  distance  of  l*  =  0(fcf 1  A) 
downstream  from  the  leading  edge.  Accordingly,  we  introduce  the  Reynolds  number  based  on  /*: 

(2.17)  R=^L=Ra{1-). 

The  base  flow  which  has  been  distorted  by  this  low-frequency  fluctuation  is  inhomogeneous  in  both 
space  and  time.  Although  no  comprehensive  theory  exists  to  describe  the  complete  instability 
characteristics  of  such  flows,  we  will  exploit  the  weakness  of  the  inhomogeneity  to  identify  some  of 
the  interesting  features  of  this  instability. 

Similar  to  the  case  of  purely  steady  distortions  due  to  streamwise  vortices  (Goldstein  & 
Wundrow  [14];  Wu  &  Luo  [40]),  the  spanwise  velocity  perturbation  induced  by  the  Klebanoff 
fluctuation  is  small  compared  with  the  corresponding  perturbation  in  the  streamwise  velocity. 
Therefore,  to  the  leading  order  of  approximation,  the  modified  base  flow  corresponds  to  a  unidi¬ 
rectional  shear  flow,  with  the  primary  shear  being  along  the  wall-normal  direction  and  a  secondary 
shear  along  the  spanwise  coordinate.  Such  flows  can  support  viscous-inviscid  interactive  instabili¬ 
ties  with  0(R1/AUOQir)  frequencies,  which  are  governed  by  the  triple-deck  structure  (Smith  [30]). 
However,  it  is  easily  verified  that  the  amplitude  of  the  Klebanoff  distortion  must  be  comparable  to 
the  free-stream  speed  in  order  to  produce  an  0(1)  change  in  the  growth  rates  of  these  instability 
modes.  In  contrast,  even  relatively  weak  Klebanoff  distortions  can  have  a  significant  effect  on  the 
instability  modes  at  higher  frequencies  (compared  with  those  of  the  viscous-inviscid  interactive 
modes),  whose  growth  characteristics  are  controlled  by  the  curvature  of  the  perturbed  flow  in  the 
vicinity  of  the  wall.  An  important  observation  in  this  context  is  that,  for  the  low-frequency  (but 
unsteady)  distortion, 

(2.18)  U"  ~  rj  as  rj  0  . 
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Therefore,  in  a  thin  region  near  the  surface,  the  Klebanoff  fluctuation  at  suitable  values  of  k\  and 
to  may  alter  the  curvature  of  the  unperturbed  Blasius  profile  by  0(1),  while  the  perturbation  to 
the  streamwise  velocity  itself  remains  small  compared  with  the  free-stream  speed.  This,  in  turn, 
can  lead  to  fundamental  changes  in  the  nature  of  instability  in  the  distorted  flow,  by  allowing 
inviscid  instability  modes  to  exist  in  addition  to  the  primarily  viscous  T-S  modes  supported  by 
the  Blasius  flow.  The  above  scenario  is  rather  different  from  the  case  of  a  completely  steady 
distortion,  for  which  U"  ~  rf  as  rj  0  and,  therefore,  no  new  instability  can  emerge  until  the 
distortion  amplitude  becomes  0(1).  Thus,  there  exists  a  crucial  difference  between  steady  and 
unsteady  distortions,  even  if  the  frequency  of  the  distortion  is  significantly  smaller  than  that  of 
the  instability  modes. 

The  curvature  of  the  distorted  streamwise  profile  becomes  comparable  with  that  of  the  Blasius 
profile  itself  within  a  thin  region  near  the  wall.  The  0(a)  thickness  of  this  region  (relative  to  the 
0(R—  i /*)  thickness  of  the  mean  boundary  layer)  is  determined  by  the  balance 


(2.19)  cjDd(^-)  ~  a2  . 

Given  the  possible  occurrence  of  an  inflection  point  at  distances  of  0(a)  from  the  surface,  it  follows 
from  the  classical  Rayleigh  scalings  that  the  most  unstable  inviscid  modes  will  have  phase  speeds  of 
0(aUoo)  and  wavelengths  of  0(d_1i?_^/*)  along  the  streamwise  direction.  Furthermore,  consistent 
with  our  focus  on  local  distortions,  we  seek  instability  modes  whose  streamwise  wavelengths  are 
shorter  compared  with  the  spanwise  length  scale  of  A. 

The  spanwise  mode  shape  of  the  above  instabilities  (and  their  growth  rate)  is  determined  by 
the  spanwise  distribution  of  the  Klebanoff  distortion  at  the  location  of  interest,  similar  to  the 
instability  modes  considered  in  Part  I  (which  involved  purely  stationary  distortions).  A  crucial 
difference  from  Part  I  is  that  the  curvature  of  the  Klebanoff  distortion  is  smaller  by  a  factor  a  in  the 
wall  layer  than  in  the  main  part  of  the  boundary  layer.  Unlike  in  Part  I,  therefore,  the  instability 
characteristics  of  the  perturbed  flow  are  not  solely  controlled  by  the  distortion  in  the  wall  region. 
The  distortion  in  the  bulk  of  the  boundary  layer  also  affects  the  instability  wave.  We  will  show  in 
§3  that  balancing  the  ‘scattering5  effect  due  to  the  spanwise  modulation  with  the  amplification  of 
of  the  instability  mode  yields  the  constraint 


(2.20) 


„  3 

a 2  ~ 


R~il* 

A 


It  is  easy  to  deduce  from  (2.17),  (2.19)  and  (2.20)  that  the  width  of  the  wall  layer  a  is  related 
to  the  Reynolds  number  R\  and  the  frequency  k\  of  the  Klebanoff  fluctuation  via 


(2.21) 


o  ~  {Rh/ki)  s  , 


and  that  the  required  amplitude  of  the  free-stream  disturbance  is 
(2.22)  to  ~  k\a  ~  Ra  3  Aq3  . 


We  thus  set 

eD  =  R^kfB0  , 


li 


so  that  Bo  =  0(1)  denotes  the  scaled  amplitude  of  the  distortion. 

It  can  be  shown  that  for  the  Klebanoff  distortion  to  induce  an  0(1)  (or  larger)  change  in  the 
viscous  growth  rate  of  the  instability  modes  of  interest  (relative  to  that  in  an  unperturbed  Blasius 
flow),  we  must  have  a  >  R~^ .  On  the  other  hand,  recall  that  the  foregoing  analysis  was  based 
on  the  assumption  that  d  «  1.  These  considerations  impose  the  following  restriction  on  the  ratio 
ki  of  the  spanwise  and  streamwise  length  scales  of  the  Klebanoff  distortion, 

(2.23)  R^1  <C  fci  <C  R~K™  . 

or,  equivalently,  on  the  range  of  streamwise  locations  where  the  analysis  in  §§3-4  is  formally  valid, 

(2.24)  ^  <C  i?A  • 

Correspondingly,  the  amplitude  of  the  free-stream  disturbance  can  be  anywhere  within  the  range 

i  _12 

(2.25)  i?A  «^<i?A23- 

While  the  Klebanoff  distortion  modulates  the  Blasius  flow  on  the  slow  variables  x  and  t  =  kit, 
the  instability  wave  oscillates  on  the  much  faster  variables  aRix  and  a2R^i.  We  thus  introduce 

(2.26)  (  =  (ax  —  aut)  , 

to  describe  the  oscillation  of  the  carrier  wave,  where  the  scaled  wavenumber  and  frequency  a  and 
u)  are  both  0(1),  and  a  expands  as 

a  =  ao  +  aai  +  +  a3as  +  ...  . 

It  turns  out  that  ao,  a±  and  a2  are  purely  real  quantities,  and  as  is  the  first  term  to  have  a 
nonzero  imaginary  part.  In  other  words,  the  growth  rate  of  these  instability  modes  is  of  0(d4), 
which  exceeds  the  viscous  growth  rate  of  the  lower  branch  T-S  modes  when  a  »  i?-1/32.  To 
facilitate  the  subsequent  study  of  nonlinear  development  of  the  above  modes,  we  choose  to  absorb 
both  amplitude  and  phase  variations  at  the  level  of  as  into  an  amplitude  function  that  varies  on 
the  scale 

(2.27)  X  =  a^Rhx  =  0(1)  . 

Since  X  is  also  much  faster  than  x  according  to  (2.21)  and  (2.23),  the  space  and  time  modulation 
of  the  distortion  can  be  treated  as  parametric  when  the  stability  of  the  perturbed  flow  is  studied. 

3.  Linear  instability  analysis.  As  in  Part  I,  the  linear  instability  problem  is  governed  by  a 
five-zoned  asymptotic  structure.  In  the  present  problem,  however,  the  Klebanoff  distortion  in  the 
main  deck  plays  a  crucial  role  in  determining  the  spanwise  mode  shape.  In  contrast,  the  relatively 
strong  distortion  in  the  wall  region  alone  was  found  to  control  both  the  mode  shape  and  the  growth 
rate  for  the  case  analyzed  in  Part  I. 
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The  appropriate  transverse  variable  in  the  main  part  of  the  boundary  layer  corresponds  to 

.  A  Ri 

y  =  —pry  =  (2x)2v  • 

The  total  streamwise  velocity  of  the  base  flow  is  then  given  by 

(3.1)  Ub  +  o-Un(y->  Z;x,i) 


where  Ud  denotes  the  normalised  distortion  profile.  For  the  most  part,  we  focus  on  the  the  case 
of  a  time  harmonic  free-stream  disturbance,  for  which  Ub  —  —w's(Z)(u e_i* +c.c.^.  However, 
the  linear  instability  analysis  presented  in  this  section  is  also  applicable  to  broadband  free-stream 
disturbances.  The  eigensolutions  in  the  main  deck  assume  the  usual  form 

(3.2)  u  =  e^A(X)^(Z)u0(y)  +  aui  +  a2u2  +  dsu3  +  . . .  j  elC  +c.c., 

(3.3)  v  =  ta^A(X)$(Z)v0(y)  +  av i  +  u2v2  +  dsv3  +  . . .  j  e1^  +c.c.  , 

(3.4)  w  =  (Z)w0(y)  +  awi  +  a2w2  +  asw3  +  . .  .  j  e1^  +c.c., 

(3.5)  p  =  ea{A(X)$(Z)p0(y)  +  dpi  +  a2p2  +  asp3  +  . .  .}  elC  +c.c.  , 

where  e  represents  the  amplitude  of  the  streamwise  velocity  fluctuation  associated  with  the  eigen- 
modes,  and  the  (normalized)  amplitude  function  A(X)  illustrates  the  growth  (or  decay)  of  the 
eigenmode  in  the  streamwise  direction. 

The  leading-order  terms  in  the  expansions  (3.2)-(3.5)  have  the  familiar  solution 

(3.6)  u0  =  U'B  :  v0  =  -i a0UB  ,  Po  =  Po  •  =  -(ia0)_1Pobr^1 

where  Po  is  a  constant.  The  second-order  terms  satisfy 

ia0ui  +v[  =  -iai  UBA<5> 

(3.7)  ia0UBui  +  UBvi  =  -la0P0A<5>  -  l(aiUB  -  (jo)UbA§  -  Rs  >  , 
c%U%A*  =  -p[ 

where  primes  denote  differentiation  with  respect  to  y  and  the  forcing  term 

(3.8)  Rs  =  (i aoUoUD+voU'D)AQ  =  i a0(U'BUD  -  UBU'D)A$ 


arises  from  the  interaction  between  the  instability  wave  and  the  distortion.  The  solution  to  (3.7) 
is  found  to  be 


(3.9) 

(3.10) 

(3.11) 


—  \  oioUd  (A<f> 


Vi  =  —  iolo^iUb  +  |i^  +  i oloPoUb  J 

= *,vB + f  f + + ^l-4® 


d  y 


}• 


Pi  =  Pi-  a, 


rv 

U* 

Jo 


UBd  y  , 
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where  3>i  and  P±  are  unknown  functions  of  Z  that  represent  the  complementary  solution  to  the 
system  (3.7).  The  inner  limits  (i.e.  y  — »  0)  of  (3.9)  and  (3.10)  are  given  by 

(3.12)  Vi  — >  i(cj  — -^-Pq)A<&  —  jiao^iA  +  i^o(AJD  —  AJoPo)^^  5 

(3.13)  fii  ^$iA  +  (-— \  +  \d-\J0Po)A$, 

a0 

respectively,  where  A  =  0.33206  x~i  denotes  the  local  wall  shear  of  the  Blasius  flow;  and 

(3.14)  A  n(Z]x7t)  =  UfD(07Z]x7t)  =  (2x)~iw's(Z)(u'(0;x^t)e~1^  +c.c.^ 

represents  the  instantaneous  wall-shear  perturbation  due  to  the  Klebanoff  fluctuation.  The  con¬ 
stant  Jo  is  defined  by 

Jo  =  ~Io  ay)^+w  ’ 

where  a  is  an  arbitrary  constant. 

In  the  upper  deck  region  corresponding  to  y  =  ay  =  0(1),  it  is  adequate  to  consider  the 
solution  for  the  perturbation  in  pressure  and  vertical  velocity.  These  expand  as 

(3.15)  p  =  ea^A(X)<f>(Z)p0  +  api  +  a2p2  +  asp3  +  . . .  j  elC  +c.c.  , 

(3.16)  v  =  ta$yA(X)$(Z)vo  +  av i  +  a2v 2  +  a3v3  +  . .  .|  e1^  +c.c.  . 

The  leading-order  pressure  perturbation  p0  and  the  first-order  correction  pi  are  governed  by 
Po,yy  -  ®oPo  =  o  ,  pim  -  alp 1  =  (2a0a1^  -  $Zz)Ap0  , 


which  yield  the  solutions 

(3.17)  po  =  Po  e-“oS  ,  Pi  =  Pi(X,  Z)  e-“»®  +A(-a^  +  ^-$zz)P0y  e~a°v  . 

Zao 

The  corresponding  solution  for  the  vertical  velocity  fluctuation  can  be  easily  derived  from  the 
vertical  momentum  equation, 


iaot’o  =  ~Po  ,  i«o^i  +  i(«i  -  w)A$vq  =  -pi  , 
which  implies  the  inner  behaviour 

(3.18)  vo—t—iPo-f  —  iPi  +  A{—  i — $  +  2$zz)-Po  as  J  ~ ^  0. 

a0  2o(q 

By  matching  the  upper-deck  solutions  for  pressure  and  vertical  velocity  with  the  outer  expansions 
of  the  main-deck  solutions  ((3.9)  and  (3.11)),  we  obtain 

(3.19)  Pi  =  Pi-  a%I2A$  , 

(3.20)  —  \Pi  +  (—  i  — +  2  $zz)Po  =  —  i  Q!oj4$i  +  (i  u>  +  i  do  J<xPo)A.<§>  , 
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where 

roo  roc  ^  \ 

h  =  l  uiiv,  ;„  =  |  (cjj-1)"- 

The  main-deck  expansions  (3.2)-(3.5)  become  nonuniform  as  y  0  and  one  must,  therefore, 
consider  the  Tollmien  layer  defined  by 


Y  =  y/a  =  0(1)  , 

where  the  base  flow  can  be  approximated  by  the  McLaurin  expansion  of  (3.1), 

(3.21)  <r(A  +  &XD)Y+a4(~Yi  +  U% (0, Z; x, f)F3)  +...  . 

Note  that  the  velocity  of  the  Blasius  flow  is  changed  only  slightly  by  the  distortion,  but  its  curvature 
is  completely  altered. 

The  solution  for  the  instability  wave  expands  as 


(3.22) 

u  =  e{A(X)$(Z)Uo  +  d-Ui  +  a2U2  +  <73P3  +  •  •  • 

|  elC  +c.c.  , 

(3.23) 

v  =  ea2{A(X)$(Z)V0  +  +  a2V2  +  a3V3  +  . 

. .  j  e1^  +c.c. 

(3.24) 

p  =  ea{A(X)<f>(Z)P0  +  &P1  +  a2P2  +  a3P3  +  . . 

.  j  e1^  +c.c. 

The  leading-order  solution  is 

Uq  —  A  ,  Vo  —  —  iaoAY  , 

and  matching  with  the  main  deck  solution  gives  Pq  =  (Xq:  as  well  as  the  leading-order  dispersion 
relation 

(3.25)  a0  =  (Aw)i. 

The  second-order  terms  in  the  Tollmien  layer  are  governed  by 
iotoUi  +  V[  =  —  iaqA-A^  , 

i^O'AK  —  cu)Ui  +  XVi  =  —  i otoi\  —  i  ol\PqA^  —  i  Qq  A2^4$T  . 

The  matching  requirement  with  the  small-?/  asymptote  of  the  main-deck  solution  (3.12)-(3.13) 
suggests  that  U\  and  V\  must  have  the  solution 

U1=  ^X  +  (-^X  +  XD-XJ0P0)A^ 

'  Vx=  {-ia0$iA  +  ia0(A13-AJoPo)A$}y 

The  above  solution  satisfies  the  continuity  equation  and  substituting  it  into  the  momentum  equa¬ 
tion  yields  the  constraint 

(3.27)  wA{$i  +  +  X  _  JoP°)j4$}  =  a°Pl  +  • 
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It  follows  from  the  relations  (3.19),  (3.20),  (3.25)  and  (3.27)  that  the  spanwise  mode  shape  $(Z) 
is  governed  by  the  Schrodinger  equation 

(3.28)  x,  t)  —  $  =  0  , 

where  we  have  set 

(3.29)  =  2XD(pa^  =  -J^%{u'e-'i+c.c.)w's(Z)  =  -i>(x,t)w's(Z), 

A  (2x)  2  A 

and 

=  4a0  jaq  -  ^ao(Joo  -  J0  -  J2)  -  u }  . 

Observe  that  as  represents  the  eigenvalue  of  the  Schrodinger  operator  with  a  potential  function  (or 
“scatter”  in  the  general  sense)  corresponding  to  the  distortion  A £>(Z;x,t)  in  the  wall  skin  friction. 
Even  though  the  interaction  between  the  instability  wave  and  the  distortion  takes  place  in  the  bulk 
of  the  flow,  the  net  effect  of  this  interaction  is  fully  characterized  by  a  local  quality  A d  with  the 
detailed  profile  of  the  distortion  being  largely  irrelevant. 

It  now  transpires  that  the  reason  for  choosing  (2.20)  was  to  ensure  a  balance  between  the 
spanwise  variation  $zz  and  the  wavenumber  correction  as$  in  (3.28).  Without  retaining 
would  be  parametrically  dependent  on  the  spanwise  variable  Z.  Any  higher  order  terms  involving 
first  and  second  order  derivatives  with  respect  to  Z  would  lead  to  secular  terms  proportional  to  x 
and  x 2,  respectively,  which  would  have  invalidated  the  entire  perturbation  scheme.  An  eigenvalue 
problem  analogous  to  (3.28)  was  previously  obtained  by  Timoshin  &  Smith  [32]  in  the  context  of 
purely  inviscid  singular  modes  induced  by  a  stationary  distortion. 

Thus  far,  we  have  not  made  any  distinction  between  local  and  periodic  distortions.  In  the 
rest  of  this  section  and  in  the  following  section,  we  shall  assume  that  A d(Z)  is  localised.  However, 
the  main  results  can  be  applied  to  the  periodic  distortion  after  rather  minor  modifications,  as 
discussed  in  §5  below.  The  Schrodinger  operator  is  well  studied  and  there  is  extensive  literature  on 
the  structure  of  its  spectrum.  For  our  purpose,  it  suffices  to  mention  that  for  a  localised  potential, 
the  Schrodinger  operator  has  a  discrete  spectrum  (if  the  potential  is  not  negative  definite),  and  the 
eigenfunctions  <£  are  real  valued  and  decay  exponentially  as  Z  ±oo.  In  this  paper,  we  normalize 
the  eigenfunction  such  that 


There  also  exists  a  continuous  spectrum  such  that  remains  finite  and  oscillatory  at  Z  =boo. 
However,  only  the  discrete  spectrum  will  be  considered  in  this  paper,  with  the  exception  of  a  few 
additional  comments  on  the  possible  relevance  of  the  continuum  modes. 

The  analysis  can  be  carried  to  higher  orders  in  a  routine  manner  (cf.  Wu,  Stewart  &  Cowley 
[41]).  Consideration  of  the  third  terms  in  the  expansion  for  each  deck  determines  a2,  but  the 
details  of  this  calculation,  as  well  as  a2  itself,  are  of  little  relevance  here  and  hence  are  omitted. 
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The  crucial  equation,  which  determines  the  modal  growth  rate,  is  obtained  by  considering  the 
fourth  term  in  the  expansion  for  each  deck.  The  final  result  is 

(3.30)  ~^2,zz  =  J_{2A^)qo  -  as}$2  -  Ax*  +  (70  +  i(Z))a*  +  i x(Z,X)  , 

where  the  function  $2  (X,Z)  represents  the  complementary  part  of  the  solution  for  U2  and  V2, 
analogous  to  $1  in  (3.26).  The  constant  70  is  defined  by 

_  7 tCq  A2 

4A  2i?^d5(2aoCo)^ 

while  7 (Z;x,t)  corresponds  to 

(3.31) 7 (Z\x,t)  =  Z-,x,t)(2x)~i  |C"/(0;x)e“ii  +c.c.}w's(Z)  =  -7  (x,t)w's(Z)  . 

The  function  x  is  real- valued  and  therefore  does  not  affect  the  leading-order  growth  rate. 

Equation  (3.30)  is  an  inhomogeneous  Schrodinger  equation.  The  standard  procedure  of  im¬ 
posing  the  solvability  condition  yields 


Ax  =  (70  +  Kd)A  , 


where 

(3.32) 


Kd  = 


dZ  . 


The  total  growth  rate  therefore  corresponds  to  (70  +  Kd),  with  Ka  being  the  ‘excess  growth  rate’ 
induced  by  the  distortion.  When  a  »  0(R~^),  the  excess  growth  rate  becomes  much  larger  than 
the  second  term  in  70,  which  corresponds  to  the  viscous  contribution  to  the  growth  rate  (Goldstein 
&  Durbin  [12]).  In  other  words,  the  instability  modes  of  interest  now  becomes  predominantly 
inviscid,  with  a  growth  rate  given  by 

(3.33)  K  =  -^  +  ^°°7(Z)$2d^. 

As  the  distortion  amplitude  is  increased  further,  the  inviscid  growth  rate  continues  to  increase 
and  for  a  >  0(,R-1/32)  it  exceeds  the  growth  rate  of  the  longer  wavelength  lower-branch  modes 
(which  are  described  by  the  triple-deck  structure  and  correspond  to  the  most  unstable  modes  of 
the  unperturbed  flow).  Clearly,  the  nature  of  the  instability  has  been  fundamentally  altered  at 
this  stage. 

In  summary,  we  have  seen  that  the  asymptotic  regime  studied  above  describes  a  continuous 
transition  as  the  distortion  amplitude  is  varied,  from  a  modified  form  of  the  short- wavelength, 
viscous  modes  in  an  unperturbed  Blasius  flow,  to  primarily  inviscid  modes  that  eventually  dominate 
the  overall  instability  of  the  perturbed  flow.  Because  the  structure  of  these  modes  can  be  localised 
in  the  spanwise  direction  and,  in  general,  is  completely  dictated  by  the  shape  of  the  Klebanoff 
distortion,  we  will  refer  to  these  modes  as  ‘localised  T-S  modes’. 
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1.2 


Fig.  3.1.  The  shape  of  the  distortion.  Curves  (l)-(3)  correspond  to  d  —  2,  4,  6 

It  may  be  seen  from  (3.28)  and  (3.33)  that  the  spanwise  shape  and  the  growth  rate  of  the 
instability  mode  are  controlled,  respectively,  by  the  wall  shear  U'D(0:  Z;  x,  t)  and  the  torsion 
£7^(0,  Z;x,t)  of  the  distortion  in  the  streamwise  velocity  profile.  For  a  time  harmonic  gust,  the 
latter  is  given  by 

(3.34)  U'jS(0,Z-,x,t)  =  (2$)-i{c7"'(0,Z;$)e-if+c.c.}<(Z)  , 

which,  in  view  of  (2.16),  is  related  to  the  wall-shear  distortion  [7^(0,  Z;x,  t)  via 

U'£(P,Z]x,i)  =  U'D(P,Z]x,i+^)  . 

However,  as  pointed  out  before,  the  original  equations  (3.28)  and  (3.32)  governing  the  instability 
modes  are  valid  for  any  arbitrary  form  of  time  dependence  for  the  free-stream  gust  as  long  as  the 
asymptotic  scalings  outlined  in  §2.2  are  satisfied. 

Watmuff  [35]  had  observed  that  when  the  amplitude  of  the  free-stream  disturbance  was  re¬ 
duced,  transition  actually  occurred  earlier.  This  is  indeed  an  anomalous  behaviour  if  the  level  of 
the  free-stream  disturbance  is  regarded  as  the  relevant  control  parameter.  However,  Watmuff  also 
pointed  out  that  the  decrease  in  distortion  amplitude  was  accompanied  by  the  location  of  its  maxi¬ 
mum  moving  closer  to  the  wall;  so  it  is  quite  possible  that  the  wall-shear  perturbation  UfD( 0,  Z\  x,  t) 
actually  increased  in  magnitude  during  that  process.  If  that  was  true,  then  the  seemingly  anoma¬ 
lous  observation  would  be  completely  ‘normal’  according  to  the  present  theory.  Unfortunately,  the 
wall-shear  perturbation  was  not  measured  in  Watmuff ’s  experiment  and  is  cumbersome  to  measure, 
in  general.  Therefore,  the  above  conjecture  cannot  be  easily  verified. 

The  spectrum  of  Schrodinger  operator  is  well  understood;  but,  in  order  to  aid  our  subsequent 
discussions,  we  first  present  solutions  for  a  Klebanoff  distortion  with  the  spanwise  distribution 

(3-35)  B(20  =  ^  , 
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where  d  is  a  constant.  A  localised  distribution  of  this  type  is  believed  to  be  appropriate  for  Klebanoff 
distortions  with  a  finite  correlation  distance  in  Z,  but  periodic  distributions  (corresponding  to  large 
coherence  in  Z)  can  also  be  easily  analyzed  as  discussed  in  §5.1  below.  Substituting  for  B(Z)  into 
(2.8)  and  evaluating  the  principal-value  integral,  we  find  that 


(3.36) 


B0(l-3Z2/d2) 

(zya2  + 1)3 


=  B0S(Z)  , 


where  we  have  set  B0  =  2'kvdo.  The  shape  function  S(Z)  is  depicted  in  figure  3.1  for  several  values 
of  d.  The  requirement  that  the  spanwise  velocity  of  the  distortion  vanishes  at  ±oo  dictates  that 
S(Z)  must  change  its  sign  in  the  spanwise  direction.  All  of  the  calculations  presented  below  pertain 
to  the  (arbitrary)  choice  of  d  =  4. 

Just  to  illustrate  the  general  behaviour  of  the  spectrum,  we  first  set  ip  =  1  in  (3.29)  and  plot 
the  eigenvalues  as  for  a  range  of  B0  (figure  3.2).  The  variation  of  B0  also  emulates  various  instances 
during  a  single  period  of  the  time-harmonic  Klebanoff  fluctuation.  Specifically,  Bo  corresponds  to 
the  modulation  phase  in  which  the  distortion  is  characterized  by  a  low-speed  streak.  Symmetric 
(varicose)  modes  can  be  found  for  both  positive  and  negative  B0:  except  in  a  gap  0  <  B0  <  0.16  = 
Bc.  Unlike  the  Schrodinger  operator  with  a  purely  imaginary  potential  (i.e.  the  case  analyzed  in 
Part  I),  the  standard  Schrodinger  equation  (3.28)  also  admits  antisymmetric  (i.e.  sinuous)  modes 
in  addition  to  the  varicose  ones.  These  modes  appear  only  for  negative  B0  (or,  equivalently,  during 
a  small  window  that  does  not  exceed  one  half  of  the  modulation  cycle  in  terms  of  duration).  For 
Bq  <  0,  there  are  also  higher  unstable  modes,  both  symmetric  and  antisymmetric,  which  are 
distinguished  by  the  number  of  zeros  in  the  corresponding  eigenfunctions,  as  displayed  in  figure 
3.3(a)  and  (b).  These  higher  modes  are  generally  less  unstable  than  the  first  ones  and,  accordingly, 
will  not  be  discussed  here  any  further.  While  the  symmetric  modes  for  B0  >  Bc  concentrate  near 
the  centerline  (see  the  dashed  line  in  figure  3.3b),  the  modes  for  Bo  <  0  have  peaks  outside  the 
main  region  of  the  mean-flow  distortion.  Interestingly  enough,  both  the  eigenvalue  as  and  the 
eigenfunctions  4>(Z)  for  the  varicose  and  sinuous  modes  for  B0  <  0  are  remarkably  similar  to  each 
other  (except  in  the  immediate  vicinity  of  Z  =  0)  in  this  particular  case,  as  seen  from  figures  3.2 
and  3.3. 

To  assess  the  relative  importance  of  the  sinuous  and  the  varicose  modes,  we  now  examine  the 
‘unit  growth  rate’,  which  is  defined  as 


«o  =  - 


S(Z)$2dZ 


and,  therefore,  represents  the  excess  growth  rate  (3.32)  normalized  by  the  amplitude  of  the  distor¬ 
tion.  As  is  shown  in  figure  3.4,  the  first  sinuous  modes  exhibit  a  slightly  higher  unit  growth  rate 
than  the  varicose  modes,  except  as  B0  — >  0  when  the  growth  rate  kq  decreases  sharply. 

Equations  (3.28),  (3.29)  and  (3.31)  are  used  in  conjunction  with  (3.33)  to  compute  the  inviscid 
growth  rate  due  to  the  time- varying  Klebanoff  distortion.  The  growth  rates  of  both  the  sinuous 
and  varicose  modes  at  three  separate  instants  of  time  are  shown  in  figure  3.5  for  the  case  of  x  =  2.0 
and  B0  =  1.4.  Observe  that  the  sinuous  modes  have  considerably  larger  growth  rates  than  the 
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Fig.  3.2.  Eigenvalues  as  versus  Bq:  solid  line:  varicose  modes,  dashed  line:sinuous  modes. 


Fig.  3.3.  Eigenfunctions  of  the  first  and  second  modes  for  the  localised  Klebanoff  distortion  defined  by  (3.36): 
sinuous  mode  (a)  and  varicose  mode  (b)  for  d  —  4  and  Bo  =  —1.  For  comparison,  the  dashed  line  represents  the 
eigenfunction  of  the  varicose  mode  for  d  —  4  and  Bq  =  1 . 
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Fig.  3.4.  Variation  of  the  ‘unit  growth  rate ’  with  Bo.  Solid  line:  varicose  modes ;  dashed  line:  sinuous  modes. 

varicose  modes.  For  this  reason,  we  shall  focus  on  the  sinuous  modes  henceforth.  Indeed,  even 
in  the  experiments  of  Matsubara  &  Alfredsson  [26],  the  sinuous  modes  were  noted  to  occur  more 
frequently. 

As  described  earlier  in  the  context  of  figure  3.2,  sinuous  modes  exist  only  when  /ip(x7i)  is  neg¬ 
ative,  i.e.  at  those  instants  during  the  modulation  cycle  when  the  perturbed  flow  is  characterised 
by  a  significant  low-speed  streak.  This  finding  is  consistent  with  many  of  the  experimental  obser¬ 
vations  noted  in  the  Introduction.  Furthermore,  the  instability  occurs  only  in  that  part  of  cycle 
when  7 (x,t)  is  also  negative.  On  the  other  hand,  7  and  'ip  tend  to  0  as  both  x  0  and  x  00; 
the  former  is  implied  by  (2.15),  while  the  latter  is  apparent  from  figure  2.1.  Thus,  7  and  pj  have 
appreciable  magnitudes  over  only  a  restricted  window  in  the  streamwise  direction.  The  instability 
modes  under  consideration  are,  therefore,  expected  to  be  localised  in  space  as  well  as  in  time. 

The  local  and  intermittent  nature  of  the  instability  can  be  further  illustrated  by  plotting  the 
growth-rate  contours  in  the  u  —  x  plane  at  various  instants  of  time,  as  shown  in  figures  3.6(a)-(d) 
for  the  case  of  Bo  =  1.4.  At  t  =  —1.8,  a  small  ‘bubble’  of  instability  is  observed  within  the  u)  —  x 
plane,  indicating  the  incipience  of  the  instability  to  be  at  a  slightly  earlier  time.  As  time  increases, 
the  bubble  grows  in  both  spatial  and  spectral  extent,  reaching  its  maximum  at  t  &  —0.82,  after 
which  the  bubble  shrinks  and  finally  disappears,  before  re-emerging  during  the  next  cycle  of  the 
Klebanoff  fluctuation.  Of  course,  given  the  disparity  between  the  temporal  scales  of  the  Klebanoff 
fluctuation  and  the  instability  waves  the  latter  could  amplify  substantially  within  a  single  period 
of  modulation  and,  therefore,  reach  sufficiently  high  amplitudes  to  induce  a  local  breakdown.  For 
the  periodically  modulated  Klebanoff  distortion,  the  instability  occurs  during  a  fixed  phase  of  the 
cycle,  but,  in  reality,  the  modulation  is  stochastic  so  that  this  local  instability  takes  place  randomly 
as  has  been  observed  in  the  experiments  of  Kendall  [18]. 

It  is  possible  to  deduce  useful  information  about  the  spatio-temporal  behaviour  of  the  instabil¬ 
ity  from  figures  3.6(a)-(d).  The  localised  structure  of  the  instability  wave  in  both  x  and  Z  suggests 
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Fig.  3.5.  Instability  caused  by  a  Klebanoff  distortion.  The  streamwise  location  is  fixed  atx  —  2.0  and  Bo  =  1.4. 
The  figure  shows  the  growth  rates  of  varicose  (solid  line)  and  sinuous  modes  (dashed  line)  at  three  instants:  (1) 
t  —  <fi(x)  —  |7r;  (2)  t  =  <f{x)  —  |7r;  (3)  t  =  <p(x)  —  ||7r,  where  Pd  2.31  denotes  the  local  phase  of  the  wall  torsion 
U'"{  0,x). 


Fig.  3.6.  Contours  of  growth  rates  of  the  local  instability  induced  by  a  Klebanoff  mode  (Bo  —  1A).  Figures 
(a)-(d)  correspond  to  the  instants  t  —  —1.8,  -1.25,-0.82,-0.44. 


that,  in  practice,  the  unstable  modes  would  be  manifested  as  patches  of  oscillations.  Strictly  speak¬ 
ing,  the  centres  of  these  patches  can  only  be  determined  after  considering  the  complex  history  of 
all  unstable  disturbances.  However,  it  may  still  be  useful  to  model  the  centre  of  a  patch  as  being 
close  to  the  x  location  that  corresponds  to  the  growth  rate  maximum  in  figures  3.6(a)-(d).  ^From 
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t  =  —1.80  to  t  —  —0.82,  the  centre  of  the  patch  migrates  from  x  =  0.8  to  1.6,  yielding  an  apparent 
patch  ‘convection  velocity’  of  Vc  &  0.82.  Even  though  the  above  prediction  is  based  on  a  number 
of  simplifying  assumptions,  migration  of  patches  at  such  a  speed  appears  to  have  been  observed  in 
experiments  as  well.  It  may  be  stressed  that  the  estimation  of  the  patch  convection  speed  became 
possible  only  because  the  unsteady  nature  of  the  distortion  was  accounted  for  with  a  reasonable 
accuracy. 

It  is  expected  that  the  unstable  modes  might  amplify  sufficiently  over  a  single  cycle  of  modula¬ 
tion,  so  as  to  unleash  a  chain  of  nonlinear  events  as  described  in  the  following  section.  However,  if 
the  nonlinear  effects  do  not  come  into  play  within  a  single  period  (or  streamwise  wavelength)  of  the 
modulation  (i.e.  starting  from  the  genesis  of  the  instability  mode),  then  the  question  arises  as  to 
how  a  given  disturbance  entity  evolves  through  multiple  modulation  cycles  in  time  and/or  in  space. 
It  may  be  possible  to  pursue  a  multi-dimensional  Floquet  theory  (i.e.  Floquet  analysis  in  both  x 
and  t)  within  a  finite-Reynolds-number  framework  in  order  to  examine  how  the  local/instantaneous 
modes  are  linked  to  the  global  (i.e.  Floquet)  ones.  Such  a  line  of  attack  has  has  recently  been  taken 
by  Hall  [17]  on  the  analogous  problem  for  time-periodic  (but  spatially  homogeneous)  Stokes  layer. 
Hall  finds  that  the  phase  speeds  of  certain  instantaneously  growing  modes,  if  continued  into  their 
decaying  phase,  have  a  periodic  dependence  on  time,  and  thus  these  modes  correspond  to  the  high- 
Reynolds-number  version  of  the  Floquet  modes.  However,  all  such  modes  are  found  to  be  globally 
stable,  which  suggests  that  local  modes  rather  than  the  global  modes  are  the  cause  of  transition,  at 
least,  in  the  Stokes-layer  case.  The  issue  of  local  versus  Floquet  modes  becomes  more  complicated 
for  the  doubly  inhomogeneous  basic  state  examined  herein.  There  is  insufficient  information  from 
the  experiments  in  terms  of  whether  transition  does  result  from  an  accumulated  growth  associated 
with  the  global  modes. 

4.  Nonlinear  instability.  The  linear  stability  analysis  in  the  previous  section  was  based  on 
the  assumption  that  the  amplitude  e  of  the  instability  waves  is  sufficiently  small.  However,  as  the 
waves  are  amplified  and  attain  a  certain  threshold  range  of  amplitudes,  the  nonlinear  effects  also 
become  significant.  In  this  section,  we  will  show  that  the  instability  waves  evolve  through  two 
distinct  weakly  nonlinear  stages  before  finally  entering  the  strongly  nonlinear  regime  of  develop¬ 
ment.  To  a  great  extent,  the  details  of  this  nonlinear  analysis  are  similar  to  the  previous  analyses 
of  nonequilibrium  critical  layers;  therefore,  many  of  the  details  related  to  this  analysis  are  omitted 
from  this  paper. 

4.1.  Nonlinear  stage  I.  The  first  nonlinear  stage  is  attained  when 
(4.1)  e  =  0(a^)~(^)“¥. 

At  this  stage,  the  fluctuations  associated  with  the  instability  wave  still  expand  according  to  (3.2)- 
(3.5),  as  in  the  linear  case.  Moreover,  the  first  three  terms  in  the  expansions  are  unaffected  by 
nonlinearity  so  that  the  shape  of  the  spanwise  eigenfunction  4>(Z)  continues  to  be  governed  by 
(3.28).  It  is  now  well  established  that  the  dominant  nonlinear  interaction  occurs  within  the  critical 
layer  and  the  effect  of  this  cubic  nonlinearity  is  to  produce  a  velocity  jump  N  across  this  layer 
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(Goldstein  [10],  Cowley  &  Wu  [7]).  The  disturbance  amplitude  in  (4.1)  was  chosen  so  that  this 
jump  influences  the  overall  development  (i.e.  growth)  of  the  instability  mode  at  the  same  order  as 
the  linear  dynamics. 

The  form  of  the  nonlinear  term  depends  on  viscous  effects.  As  in  Part  I.  we  assume  that 

R~i  =ra 13  , 

where  the  Haberman  parameter  r  (r  =  0(1))  characterizes  the  relative  influence  of  viscosity  with 
respect  to  nonlinearity  within  the  critical  layer. 

The  jump  N  is  given  by  the  same  expression  as  (5.3)  in  Part  I,  provided  that  A(X)  in  that 
equation  is  replaced  by  the  combined  mode  shape  A(X)$(Z).  Of  course,  the  Z-dependent  factors 
can  now  be  brought  outside  of  the  integral  sign.  Again,  matching  the  solution  from  each  deck  to 
the  respective  solutions  from  the  adjacent  regions  yields  (3.30),  but  with  the  nonlinear  velocity 
jump  N  added  to  the  right-hand  side.  Applying  the  solvability  condition  to  this  inhomogeneous 
equation,  we  obtain  the  nonlinear  amplitude  evolution  equation 

rOC  rOO 

(A.2)Ax  =  kA-{T  /  e(^  +  r,)K(^r,\s)A(X-0A(X-^-r,)A*(X-2^-r,)dUv  , 

Jo  Jo 

where 

/OO  1 

d  Z  ,  K(Z,  rj\s)  =  e~s^3+3^  ,  s  =  -c^A2r. 

-OO  ^ 

The  appropriate  initial  condition  for  (4.2)  follows  from  matching  the  upstream  behaviour  of  the 
nonlinear  solution  with  the  linear  stage,  namely, 

(4.4)  A  ->>  eKX  as  x  —  oc  . 


Numerical  solutions  presented  below  demonstrate  that,  for  relatively  weak  viscosity  (i.e.  r  <C 
1),  A(X)  develops  a  singularity  at  a  finite  location  Xs,  and  that  the  structure  of  this  singularity 
is  the  same  as  that  proposed  earlier  by  Goldstein  &  Choi  [11]: 


a0 

(Xs-X)3+i* 


X^XS1 


where  q  is  a  real  number.  The  above  singularity,  however,  does  not  occur  when  the  viscous 
parameter  s  exceeds  a  critical  value. 

Equation  (4.2)  is  formally  derived  by  assuming  that  r  =  0(1).  However,  the  permissible  range 
of  r  is  in  fact  quite  large, 


R 


-l 

A 


<r<  R^3  , 


as  may  be  deduced  from  (2.17),  (2.21)  and  (2.23).  Thus,  r  (or,  equivalently,  s)  can  be  either  a  small 
or  a  large  parameter  without  invalidating  (4.2).  The  simplified  form  of  this  evolution  equation  in 
the  very  viscous  limit  s  oc  can  be  derived  by  rescaling  the  amplitude  function  according  to 

(4.6)  A  =  ssi(X)eisie  . 
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Inserting  this  into  the  nonlinear  term  in  (4.2),  and  performing  integration  by  parts,  we  obtain 

r oc  poo 

N  =  (6s)-1  J  j  K(S, »?){l-  sift# (X-0  +  Q'(X-Z-r1)-2Q\X-2Z-r1)}} 
exp{iS3[0(x  -  0  +  e(X  -  Z  -  ri)  +  Q(X  -  2£  - 1})]} 

(4.7)  A(X-t)A(X-t-tl)A*(X-2£-rl)dZdTi  +  -~  * 

where  we  have  ignored  the  terms  which  do  not  affect  the  first  two  orders  of  the  asymptote  for  N. 
As  in  Wu,  Stewart  &  Cowley  [41],  we  now  introduce  the  substitution  £  — >  and  take  the  limit 

s  — y  oo  to  obtain  a  two-term  asymptotic  approximation  for  N.  Inserting  this,  along  with  (4.6), 
into  (4.2),  we  find  that 

0'  =  jfVM(*  -  v)\2 dr,  , 

(4.8)  A'  =  kA+-A  J '  ri-1[&(X)-&(X-ri)\\A(X-rj)\2dri- 

The  first  of  these  equations  indicates  that  nonlinearity  induces  a  wavelength  shortening  or  ‘dilation,5 
which  in  turn  influences  the  energy  growth  of  the  disturbance  via  the  coupling  between  the  phase 
and  the  amplitude  equations.  The  energy  growth  is  now  governed  by  an  equation  with  a  quintic 
nonlinearity,  as  may  be  verified  by  substituting  the  first  equation  into  the  second.  An  analogous 
finding  was  noted  earlier  by  Wu,  Leib  &  Goldstein  [42]  in  the  context  of  nonlinear  interaction 
between  pairs  of  T-S  waves. 

4.2.  Nonlinear  stage  II.  The  structure  (4.5)  of  the  singularity  during  stage  I  of  the  nonlinear 
evolution  shows  that 

Ax/A-iXs-X)-1  as 

which  suggests  that  the  instability  mode  must  now  evolve  over  a  faster  scale  in  the  streamwise 
direction.  In  particular,  a  distinct  asymptotic  regime  is  reached  when  Xs  —  X  =  0(d2),  because 
the  amplitude  growth  rate  becomes  comparable  with  the  0(a)  wavelength  correction  caused  by  the 
Klebanoff  distortion.  The  scaled  streamwise  coordinate  governing  this  second  stage  of  nonlinear 
evolution  is  defined  as 

(4.9)  X  =  (X-  Xs)/a2  . 

The  disturbance  amplitude  tA(X)  has  increased  to  0(e),  where 

(4.10)  e  =  ea~6  =  . 

In  this  regime,  the  instability- wave  perturbations  in  the  main  part  of  the  boundary  layer  expand 
according  to 

u  =  e|^A(-A,  Z)uq  e^  -\-aui  -1- . . .  +  c.c  , 
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with  similar  expansions  for  v ,  w  and  p.  Observe  that,  unlike  in  the  preceding  stages  of  evolution,  the 
growth  of  the  unstable  fluctuations  is  nonuniform  across  the  relevant  range  of  spanwise  locations; 
consequently,  the  amplitude  function  A  depends  on  both  X  and  Z. 

Matching  the  solutions  in  the  various  decks  at  the  second  order  of  approximation  leads  to  the 
amplitude  equation  for  A(X,  Z): 

(4.11)  A*  -  ^-Azz  =  -as)A  +  itN(X,Z)  , 

where  F  =  7tA3Qq£:4,  and  the  nonlinear  term, 


N 


(4.12) 


POO  pOO  ( 

J  J  U3A(X-0A(X-t-r))A*zz(X-2{-r)) 

H2vAX  -  0  [MX  -  i  -  rj)Az(X  -  2£  -  r,)}  z 
+£3  [A(x  -  oMx  - 1-  -  v)A*z(x  -  2^-»?)]z|d^d»j, 


again  represents  the  effect  of  interactions  within  the  critical  layer.  Because  of  the  shorter  streamwise 
length  scale  during  the  second  nonlinear  stage,  the  width  of  the  critical  layer  has  increased  from 
an  0(d4)  value  (during  the  first  nonlinear  stage)  to  0(a2).  As  a  result  of  the  thicker  critical-layer 
region,  viscous  effects  no  longer  play  any  role  in  determining  the  disturbance  evolution  at  the  above 
order. 

Equation  (4.11)  reflects  the  balance  between  spanwise  scattering  by  the  Klebanoff  distortion 
and  the  nonlinear  effects,  while  the  linear  growth  has  become  a  secondary  effect  during  stage  II.  A 
distinguishing  feature  of  this  stage  is  that  the  spanwise  distribution  of  the  instability  wave  is  altered 
by  the  nonlinear  effects,  in  contrast  to  stage  I  wherein  the  spanwise  modal  shape  was  independent 
of  the  modal  amplitude. 

The  appropriate  initial  condition  for  (4.11)  follows  from  the  requirement  of  matching  with  the 
singular  solution  from  the  upstream  region.  Thus, 

(4.13)  i-^a0(-^)_(3+i9){$(2)  +  (-X)-1$i(Z)  +  ...}ei“^  as  X  ->  -oo  . 


The  dominant  term  on  the  right-hand  side  of  (4.13)  is  easily  derived  by  rewriting  (4.5)  in  terms  of 
X.  The  governing  equation  for  the  second-order  correction  term  can  be  obtained  by  inserting 
(4.13)  into  (4.11)  and  equating  the  terms  at  0((—X)~4).  To  obtain  a  unique  solution  for  from 
this  equation,  the  orthogonality  condition, 


4>4>i  d  Z  =  0 


is  imposed.  Note  that  the  amplitude  evolution  equation  (4.11)  for  stage  II  could  have  been  derived 
without  considering  the  solution  from  stage  I.  However,  the  latter  is  still  essential  to  derive  the 
initial  condition  (4.13). 
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Fig.  4.1.  Nonlinear  evolution  of  the  most  unstable  (sinuous)  mode  in  figure  3.6b  (Bo  =  1.4,  x  —  1.2,  oj  =  0.65,) 
for  viscosity  parameter  5  =  0,  1.0,  2.0  and  15.0.  The  dashed  line  represents  the  local  singular  solution  (4.5). 


4.3.  Numerical  results  for  nonlinear  evolution.  To  illustrate  the  features  of  amplitude 
evolution  during  nonlinear  stage  I,  we  plot  the  magnitude  of  A  as  a  function  of  the  streamwise 
coordinate  X  (figure  4.1).  The  parameter  x,  which  indicates  the  streamwise  location  on  the  longer 
streamwise  scale,  is  chosen  (somewhat  arbitrarily)  to  be  1.2;  similarly,  the  distortion  amplitude  is 
chosen  to  be  B0  =  1.4  and  the  frequency  of  the  most  unstable,  i.e.  sinuous  mode  to  be  uj  =  0.65. 
For  the  purpose  of  plotting  the  results,  we  have  renormalized  the  amplitude  function  A(X)  such 
that  the  coefficient  T  =  1. 

As  mentioned  earlier,  the  solution  develops  a  singularity  of  the  form  (4.5)  at  sufficiently  small 
values  of  the  viscosity  parameter  s.  Figure  4.1  shows,  however,  that  the  amplification  is  not 
monotonic,  rather  it  exhibits  transient  decay  or  oscillations  before  the  eventual  blowup  at  a  finite 
distance  downstream.  It  is  interesting  to  note  that  increasing  the  viscous  effect  appears  to  produce 
stronger  oscillations  upstream  of  the  singularity.  The  amplitude  curve  for  s  =  1  displays  multiple 
spikes  even  on  the  logarithmic  scale.  These  spikes  are  adequately  resolved  during  the  numerical 
integration  process  by  using  a  fine  step  size  (AX  =  1/400).  At  larger  values  of  the  viscosity 
parameter  s  (viz.,  at  s  =  2  and  s  =  15),  nonlinearity  has  a  stabilising  effect  on  the  disturbance 
evolution  and  the  solution  either  decays  in  a  non-monotonic  fashion  (as  at  s  =  2  in  figure  4.1)  or 
amplifies  slowly  towards  the  end  of  the  integration  domain  (s  =  15). 

As  noted  earlier,  the  stage  I  formulation  becomes  invalid  when  approaching  the  singularity  and 
one  needs  to  solve  (4.11)  to  predict  the  subsequent  evolution  of  the  instability  mode  during  stage 
II  of  nonlinear  evolution.  In  order  to  ensure  a  smooth  match  with  upstream  solution,  however,  it 
is  necessary  to  pick  a  large  enough  value  of  —X  for  imposing  the  constraint  (4.13).  To  this  end, 
we  perform  a  trial  integration  through  the  far-upstream  region  (typically  starting  at  X  =  —250) 
so  as  to  identify  the  region  over  which  the  one-  and  two-term  approximations  from  (4.13)  may  be 
considered  as  being  sufficiently  accurate.  A  relatively  coarse  grid  (AX  =  0.5)  is  typically  adequate 
for  this  purpose.  Typically,  we  found  that  a  two-term  approximation  is  accurate  up  to  X  «  —120. 
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Fig.  4.2.  Nonlinear  evolution  in  stage:  E  v.s.  X. 


The  solution  downstream  from  this  location  is  obtained  by  marching  (4.11)  with  a 
size.  Because  the  amplitude  function  now  depends  on  both  streamwise  and  spanwise 
we  introduce  the  norm 

Uoo  1 

JA(Z,X)\2dZ}2 

to  measure  the  disturbance  amplitude  at  a  given  streamwise  location.  Figure  4.2  illustrates  the 
streamwise  development  of  E(X)  through  the  nonlinear  stage  II.  Observe  that  the  simultaneous 
spatio-temporal  modulation  during  nonlinear  stage  II  has  a  strong  destabilising  effect,  which  causes 
an  explosive  amplification  of  the  disturbance  to  lead  to  yet  another  singularity  at  a  finite  distance 
Xs  on  the  shorter  streamwise  scale  involved.  The  spanwise  distribution  of  A  at  three  typical 
streamwise  stations  is  shown  in  figure  4.3.  As  illustrated  therein,  the  nonlinear  deformation  first 
leads  to  a  shift  in  the  peak  location  from  Z  «  12.0  to  Z  &  6.5,  and  then  to  a  steepening  of  the 
mode  shape  near  the  peak  at  Zs  «  6.5.  Thus,  nonlinearity  inhibits  any  lateral  spreading  and, 
apparently,  leads  to  a  singularity  of  self-focusing  type. 

We  conjecture  that  the  structure  of  the  self-focusing  singularity  may  be  of  the  same  type  as 
that  proposed  in  Wu  [38]  for  a  similar  evolution  equation,  namely, 

(4.14)  A  =  (Xs-X)-^+ib^F(Z)  ,  Z  =  (Xs-X)~i(Z-Zs)  , 

where  b  is  real  constant  and  F  satisfies  a  rather  complex  equation  as  given  in  Wu  [38].  If  the 
final  outcome  of  the  stage  II  evolution  is  represented  by  (4.14),  then  the  subsequent  disturbance 
evolution  must  occur  on  the  even  shorter  scale 


smaller  step 
coordinates, 


-  x  =  O(d) . 


The  growth  rate  of  the  disturbance  is  now  comparable  with  the  streamwise  wavenumber  and  the 
local  spanwise  length  scale  (near  Z  =  Zs)  is  also  of  the  same  order  as  the  streamwise  wavelength. 
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Fig.  4.3.  Spanwise  distribution  of  \A\  at  different  values  of  X. 

The  appropriate  local  variables  are  then  given  by 

X  =  (X-X8)/&,  Z  =  a~HZ-Zs )  . 

One  also  finds  that  due  to  a  further  increase  in  the  thickness  of  the  nonequilibrium  critical  layer 
in  the  new  regime,  the  critical  layer  actually  merges  with  Tollmien  layer.  Since  the  amplitude 
of  the  disturbance  increases  to  —  0(d),  i.e.  is  of  the  same  order  as  the  Blasius  flow  in  the 
Tollmien-layer,  the  flow  in  this  region  becomes  strongly  nonlinear.  In  fact,  the  whole  flow  is  now 
described  by  the  inviscid  form  of  the  fully  three-dimensional  triple-deck  system.  Earlier,  Wu  et  al. 
[42]  had  shown  that  pairs  of  viscous  T-S  waves  in  the  Blasius  boundary  layer  may  evolve  through 
several  weakly  nonlinear  regimes  before  ending  up  with  the  same,  fully  nonlinear,  inviscid  triple¬ 
deck  system.  A  final  noteworthy  feature  of  the  overall  nonlinear  development  is  that  the  Klebanoff 
distortion  exerts  its  influence  only  during  the  stages  I  and  II  and  that  its  effect  is  felt  only  as  a 
higher-order  correction  during  stage  III.  In  other  words,  the  Klebanoff  distortion  is  likely  to  serve 
only  as  a  catalyst  that  promotes  the  early  growth  of  the  disturbance,  which  eventually  evolves 
into  the  same  canonical  regime  of  inviscid,  strongly  nonlinear  behaviour  as  that  predicted  for  an 
unperturbed  Blasius  flow. 
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We  emphasize  that  the  results  of  the  nonlinear  theory  ought  to  be  treated  with  some  cau¬ 
tion.  While  the  prediction  that  the  evolution  would  occur  over  progressively  shorter  length  scales 
is  broadly  in  agreement  with  experimental  observations,  the  formation  of  the  final  singularity  is 
certainly  non-physical.  Probably  strong  nonlinearity  would  prevent  it  from  occurring.  Another 
effect  that  may  alter  the  nonlinear  behaviour  is  nonparallelism  of  the  basic  state,  which  is  sup¬ 
pressed  when  integrating  the  nonlinear  evolution  equations.  Even  though  this  effect  is  small  in  the 
strict  asymptotic  sense,  it  may  be  significant  for  the  moderate  Reynolds  numbers  encountered  in 
practice. 

5.  Spanwise  periodic  distortion.  Having  examined  the  case  of  spanwise  localised  Klebanoff 
distortions  thus  far,  we  now  briefly  consider  the  instability  characteristics  in  the  presence  of  global, 
but  spanwise  periodic  distortions.  Specifically,  in  §5.1,  we  apply  Floquet  theory  to  study  the 
instability  modes  first  considered  in  §3.  In  §5.2,  we  examine  the  analogy,  for  a  nonstationary 
distortion,  of  the  Goldstein- Wundrow  [14]  (GW)  modes,  which  are  associated  with  a  subharmonic 
parametric  resonance.  The  relation  between  the  Floquet  and  the  GW  modes  is  discussed  in  §5.3. 

5.1.  Floquet  modes.  Let  us  now  examine  how  the  results  in  §§3-4  need  to  be  modified  in 
the  case  of  a  spanwise  periodic  wall-shear  perturbation  A £>(Z;x,t).  As  mentioned  earlier,  (3.28) 
still  holds  in  this  case,  but  its  solution  according  to  Floquet  theory  takes  the  form 

(5.1)  $  =  4ei/xZ, 


where  [jl  is  a  real  number,  with  /i  =  0  denoting  a  fundamental  mode  and  /i  =  1  being  the  subhar¬ 
monic  mode.  The  periodic  function  $  is  governed  by 

(5.2)  £(I>  =  &zz  +  2  i  fi&z  —  fi2&  +  , 


where  'ip(Z)  is  related  to  the  wall-shear  fluctuation  as  defined  in  (3.29).  The  operator  C 
is  Hermitian,  so  that  its  adjoint  is  simply  its  complex  conjugate,  £*.  It  also  follows  that  the 
adjoint  eigensolution  is  correspondingly  given  by  the  complex  conjugate  of  4>.  In  view  of  the 
Hermitian  property,  the  second  order  term  in  the  perturbation  expansion  (3.2)-(3.5)  merely  leads 
to  a  correction  to  the  streamwise  wavenumber  similar  to  that  for  the  localised  distortion.  At 
the  fourth  order,  again,  one  finds  (3.30)  as  obtained  earlier  in  §3.  On  writing  <f>2  =  $2 
multiplying  by  on  both  sides  of  the  resultant  equation;  and  integrating  the  resulting  equation 
over  the  spanwise  wavelength  T,  we  obtain  the  amplitude  evolution  equation  Ax  =  (70  +  Kd)A 
where  the  excess  growth  rate  Ka  is  given  by 


(5.3) 


Kd 


=  /%(Z)|$|2dZ, 
Jo 


provided  that  is  normalized  such  that  /  |4>|2  d  Z  =  1.  The  function  y(Z)  is  related  to  the  wall 

Jo 

torsion  of  the  Klebanoff  distortion  profile  (see  (3.31)). 
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Fig.  5.1.  ‘ Unit  growth  rates’  of  Floquet  modes  v.s.  fi.  Left:  branch-I  modes ;  Right:  Branch-II  modes.  Curves 
(l)-(3)  in  each  figure  correspond  to  Bq  —  0.5,  1,  2. 

Similar  to  that  in  Part  I,  a  perturbation  analysis  can  be  carried  out  for  the  limiting  case  of 
small  distortion  amplitudes  ^«1,  and  this  leads  to 

(5.4)  as  =  —l±  (<£?  +  4/x2/?2 )  i  ,  $  =  cos (3Z  +  0(ip)  , 

R  [f 

with  f3  =  'k/L  and  qc  =  —  /  i’(Z)  cos(2f3z)  d  Z.  Inserting  into  (5.3)  to  calculate  the  excess 

71  J  o 

growth  rate,  one  finds  the  total  growth  rate 

(5.5)  nr  =  Kd  +  7o  =  —  f13  'j(Z)  cos 2(/3z)  d  Z  -  ^  . 

K  JO  dA 

As  an  example,  the  calculations  will  be  carried  out  for  the  simple  harmonic  potential  function 

'ip(Z)  =  B0cos(2f3Z) 

where  Bq  denotes  a  measure  of  the  distortion  amplitude.  Unlike  the  localised  modes,  Floquet 
modes  exist  for  both  positive  and  negative  B0 ,  but  it  is  sufficient  to  just  consider  the  case  of 
Bq  >  0  since 

(5.6)  as(-B0)  =  as(B0 )  ,  4(2,  -B0)  =  4(2  +  ^,B0)  . 

The  small  amplitude  limit  (Bo  — >  0)  for  the  simple  harmonic  potential  corresponds  to 

(5.7)  as  =  —  1  =b  (-^-  +  4 ft2(32) 2  , 

which  implies  that  two  separate  branches  of  unstable  eigensolutions  exist.  Analytic  continuation  of 
these  two  branches  to  finite  values  of  B0  and  /i  will  be  referred  to  as  branches  I  and  II,  respectively. 
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The  relative  importance  of  the  two  branches  may  be  assessed  by  considering  the  ‘unit  growth 
rates’,  defined  via 

rf 

Ko=  cos(2/?Z)|$|2  dZ  • 

Jo 

The  variation  of  Ko  with  fi  is  shown  in  figures  5.1(a)-(b)  for  three  values  of  the  distortion  amplitude 
B0.  At  each  of  those  amplitudes,  the  subharmonic  mode  (/i  =  0)  has  the  largest  growth  rate,  but 
the  peak  becomes  less  prominent  as  B0  increases.  The  branch-I  modes  are  generally  more  unstable 
than  the  branch-II  modes.  Figure  5.1(b)  suggests  that  the  unit  growth  of  the  subharmonic  mode 
along  branch  II  is  a  decreasing  function  of  the  amplitude  parameter  B0. 

The  relations  (5.6)  imply  that  the  mode  for  a  negative  B0  corresponds  to  a  7t/(2/3)  shift  in 
Z  of  the  mode  for  a  positive  B0.  Since  B0  changes  its  sign  during  each  cycle  of  the  Klebanoff 
fluctuation,  the  streak  oscillation  associated  with  the  instability  is  expected  to  meander  in  the 
spanwise  direction. 

It  can  be  shown  that  the  nonlinear  amplitude  equations  (4.2)  and  (4.11),  which  were  originally 
derived  for  a  localised  distortion,  are  also  valid  for  the  case  of  a  periodic  distortion,  except  that  the 
integral  in  the  definition  of  T  (see  4.3))  now  is  over  a  single  wavelength  in  the  spanwise  direction. 

5.2.  Goldstein- Wundrow  (GW)  modes.  As  alluded  to  in  the  Introduction,  a  periodic 
distortion  may  also  be  unstable  to  the  GW  modes  (Goldstein  &  Wundrow  [14])  at  distances  Z*  = 
0(fcj“1A)  from  the  leading  edge.  Similar  to  the  local  modes  considered  in  §§3-4,  the  GW  modes 
arise  because  the  curvature  of  the  distortion  is  comparable  with  that  of  the  Blasius  profile  in  a 
wall  layer  corresponding  to  y  ~  d,  such  that 

(5.8)  cpd k^1  ~  a2  . 


The  GW  modes  also  have  0(d_1)  streamwise  wavelengths  relative  to  the  boundary-layer  thickness 
Unlike  the  Floquet  modes,  however,  both  spanwise  and  streamwise  length  scales  of  the 
GW  modes  are  comparable  with  the  spanwise  period  of  the  distortion.  Thus 


(5.9) 


R~il* 

A 


~  u  . 


The  relations  (5.8)-(5.9)  imply  that  the  threshold  disturbance  amplitude  required  for  the  onset  of 
the  GW  modes  is  given  by 


(5.10) 


CD  ~  ( 


l*R\ 

A 


{Rh/ki)  1/2  , 


Accordingly,  we  write 


ec  =  (RA/hy^Bo  , 


where  Bo  =  0(1)  denotes  the  scaled  amplitude  of  the  distortion.  A  GW  mode  consists  of  a  pair  of 
oblique  waves,  and  it  can  be  expanded  in  the  same  form  as  the  Floquet  modes,  i.e.  the  main-deck 
eigensolutions  are  given  by  (3.2)-(3.5),  except  that  in  (3.4)  is  replaced  by  d,  and 

(5.11)  a  —  ,  A$(Z)  =  eaX  (e{Pz  z be-i/3Z^  . 
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The  =b  signs  signify  two  separate  families  of  instability  modes.  Even  though  modes  from  these  two 
families  are  trivially  related  to  each  other  via  a  shift  in  Z  by  7r/(2/3),  modes  from  both  families  must 
be  considered,  in  general.  Since  a  GW  mode  involves  subharmonics  of  the  Klebanoff  distortion,  its 
amplification  can  be  attributed  to  a  parametric  resonance  mechanism. 

Similar  to  the  local  and  the  Floquet  modes,  the  Klebanoff  distortion  in  the  bulk  of  the  boundary 
layer  leads  to  an  0(a)  correction  to  the  phase  speed.  However,  for  the  present  periodic  case,  the 
above  process  is  completely  passive  in  the  sense  that  it  does  not  affect  either  the  growth  rate 
or  the  shape  of  the  instability  mode.  The  instability  characteristics  are  solely  controlled  by  the 
curvature  of  the  distortion,  exactly  as  in  Goldstein  &  Wundrow  [14].  Substituting  the  curvature 
of  the  present  distortion  into  their  equation  (6.22),  we  obtain  the  total  growth  rate  of  a  GW  mode 
induced  by  a  Klebanoff  distortion, 


(5.12) 


2  cos# 

1  +  cos2  8 


7 rex 


-AO/  -xq/o'  cos  2  0Qc(x,F)  —  H - T - r 

A2(2x)3/2  4A  2i?4d-5(2aoCo)2 


A2 


r) 


where  8  =  tan  1  — ,  c0  =  and 

cko  *  u  A  ' 

n  n^Z. 

(5.13)  nc(x,t)  =  -{^J  w's(Z)  cos(2P  Z)  d  z}  {u'"  (0,  x)  e~ +c.c.}  . 

For  the  excess  growth  rate  induced  by  the  Klebanoff  distortion  to  be  dominant,  the  disturbance 
amplitude  must  satisfy  a  >  R~^ .  Additionally,  we  must  have  a  <C  1  in  order  for  the  foregoing 
analysis  to  be  valid.  The  above  two  considerations  imply  that  the  result  (5.12)  is  valid  if  the 
wavelength  of  the  Klebanoff  fluctuation  falls  in  the  range 

(5.14)  Rj^  R\  - 

or,  alternatively,  if  the  amplitude  of  the  free-stream  disturbance  lies  within 

(5.15)  RX1  <eD  <Ra^  • 


The  result  (5.12)  indicates  that  omission  of  one  of  the  modes,  for  example,  the  minus  mode, 
will  lead  to  the  erroneous  conclusion  that  a  distortion  with  positive  >  0  will  not  be  able  to 
support  (inviscid)  instability  waves  with  cos  28  <  0.  Again,  for  a  Klebanoff  distortion,  Vtc  changes 
its  sign  during  the  modulation  cycle  and,  therefore,  there  must  be  a  switching  between  the  plus 
and  minus  modes  during  this  cycle. 

Calculations  were  performed  for  wfs(Z )  =  Bocos(2/3Z )  and  selected  values  of  (3.  In  figure  5.2, 
we  show  the  inviscid  growth  rate  (i.e.  Kr  with  the  last  term  in  (5.12)  dropped),  as  a  function  of  the 
wavenumber  ao-  For  small  and  moderate  values  of  /3,  the  peak  growth  rate  occurs  near  ao  =  0.45. 
The  peak  growth  rate  decreases  and  is  eventually  stabilised,  as  f3  is  increased.  For  larger  values  of 
/?,  the  dominant  instability  modes  emerge  at  small  a  values,  centred  at  ao  ~  0.12. 

Finally,  we  note  that  once  the  amplitude  of  the  instability  modes  reaches  the  threshold  range 
of 
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Fig.  5.2.  Growth  rates  of  Goldstein- Wundrow  modes.  The  parameters  Bq  —  A,  x  —  1.2  and  t  &  —1.498. 
Curves  (l)-(5)  correspond  to  (3  —0.063,  0.125,  0.16,  0.18,  0-44- 

the  function  A(X,  Z )  must  be  written  as 

A  =  A(X)^ei0Z±e~i0Z^  , 

where  the  amplitude  function  A  is  governed  by  the  nonlinear  equation  derived  previously  by  Gold¬ 
stein  &  Choi  [11]  and  Wu,  Lee  &  Cowley  [39]. 

5.3.  The  relation  between  Floquet  modes  and  GW  modes.  The  GW  modes  and  Flo- 
quet  modes  considered  in  §5.1  can  exist  only  when  the  Klebanoff  distortion  is  periodic  in  the 
spanwise  direction.  Let  us  discuss  the  relation  between  them  from  two  separate  viewpoints.  First, 
consider  the  stability  properties  at  a  fixed  location  of  x  =  0(1),  and  examine  which  type  of  modes 
emerge  and  how  their  character  changes  as  the  amplitude  parameter  to  is  gradually  increased  from 
zero.  As  indicated  by  (2.22)  and  (5.10),  GW  modes  require  a  smaller  threshold  for  the  distortion 
amplitude  than  the  Floquet  modes,  and  hence  the  GW  modes  appear  first  as  cd  is  increased.  That 
is,  for  to  specified  by  (5.10),  the  GW  modes  correspond  to  the  only  inviscid  instability  that  the  dis¬ 
turbed  flow  can  support.  As  to  is  further  increased,  the  streamwise  wavenumber  and  growth  rate 
of  GW  modes  (which  are  of  4)  and  O^cd/s]-1)4),  respectively)  increase  as  well.  Because  the 

spanwise  wavenumber  must  remain  fixed,  the  GW  modes  become  progressively  Two-dimensional’. 
To  investigate  this  large-amplitude  limit  (viz.  B0  >  1  or  cd  >  (Ra /&i)-1/2),  we  introduce  the 
rescaled  wavenumber  and  phase  speed  via 


(ao,  Co)  —  L?o(do,  cq)  . 


Inserting  these  into  (5.12)  and  taking  the  limit  Bo  — )>  oc  reduces  (5.12)  to 
(5-!6)  Kr  =  Bo{±a2(2^ °)3/2^($,i)  -  , 
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which  corresponds  to  the  small  distortion  amplitude  limit  of  the  growth  rate  of  Floquet  modes. 
More  precisely,  the  large-amplitude  behaviour  (5.16)  remains  valid  until  Bo  =  O^i^a)1/6),  i.e. 
cd  ~  0(R^1^3k^3),  which  corresponds  to  the  threshold  magnitude  of  the  distortion  for  the  Floquet 
modes  to  arise  (see  (2.22)).  At  this  stage,  the  enhanced  streamwise  growth  of  the  instability  modes 
is  balanced  by  the  spanwise  modulation  in  the  main  part  of  the  boundary  layer.  Thus,  we  have 
shown  that  the  GW  and  the  Floquet  modes  exist  in  separate  regimes  of  distortion  amplitude, 
and  that  the  large-distortion  limit  of  the  GW  modes  ‘matches’  to  the  small-distortion  limit  of  the 
Floquet  modes. 

An  alternative  way  to  view  the  relation  between  GW  and  Floquet  modes  is  to  consider  the 
instability  characteristics  for  a  fixed  distortion  strength  ec>,  and  examine  how  these  modes  emerge 
and  evolve  as  the  streamwise  location  is  varied.  Suppose  that  to  rsj  0(RA1^3k^3)  as  specified  by 
(2.22),  so  that  the  distorted  flow  supports  Floquet  modes  at  x  =  0(1).  Given  that  the  threshold 
distortion  strength  for  GW  modes  is  much  smaller,  one  may  expect  them  to  exist  in  the  upstream 
region  x«l.  The  link  between  the  upstream  GW  modes  and  Floquet  modes  is  best  revealed  by 
examining  the  small-x  asymptote  of  the  Floquet  modes.  If  one  takes  into  account  the  streamwise 
variation  of  the  distortion  as  well  as  of  the  Blasius  flow,  the  key  balance  (2.19)  can  be  expressed 
more  precisely  as 

(5.17)  eDk^UD,YYY(0,x)Yc  ~  a\2Y2  , 

where  Yc  =  co/A  =  ao/A2  is  the  scaled  critical  level.  Since  Ud,yyy  ~  (x)-3/2?7///,  and  U,n(0,x)  — > 
(— iAox2)  (as  indicated  by  (2.15)  and  (2.16)),  the  relation  (5.17)  shows  that  Yc  ~  x3/2  and  hence 
a o  ~  x1/2,  Co  ~  x,  as  x  0.  This  suggests  the  rescaling 

(5.18)  (ao,  co)  =  (x1/2d0,  xc0)  . 


In  the  limit  x  — >  0,  the  potential  'ip  in  (5.2)  becomes  vanishingly  small,  and  so  the  result  (5.5) 
applies.  Substituting  (5.18)  into  (5.5),  and  using  (3.29)  and  the  fact  that  U,n( 0,x)  — )>  —  iAox2,  we 
find 


(5.19) 


A223/2 


nc 


4A0  J 


where  Aq  =  0.33206,  and 


==  -  {  ^ w's(Z) cos(2 PZ)  d Z } (-  i  Ao  e- : 1  ^ "  +c.c.)  . 

The  above  approximation  becomes  invalid  as  x  =  0(fcii?A)_1//3,  at  which  stage  the  stream- 
wise  wavenumber  becomes  comparable  with  the  spanwise  wavenumber  so  that  the  modes  be¬ 
come  fully  three-dimensional  in  character.  Thus,  the  GW  modes  operate  in  the  region  where 
x  =  (kiR^Y^x  =  0(1),  and  their  growth  rates  are  given  by  (5.12),  provided  x  is  replaced  by  x, 
and  U/n  in  (5.13)  by  its  upstream  asymptote  (—  ix2Ao).  Clearly,  the  large-x  limit  of  this  result 
matches  to  (5.19).  Thus  the  upstream  limit  of  Floquet  modes  matches  to  the  downstream  limit  of 
the  GW  modes,  i.e.  the  Floquet  modes  represent  the  downstream  continuation  of  the  GW  modes. 
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One  must  note  that  the  above  matching  between  the  GW  and  the  Floquet  modes  was  based 
on  the  assumption  of  linear  instability,  at  least,  through  the  region  of  overlap  between  the  two.  In 
reality,  of  course,  the  upstream  GW  modes  might  reach  a  nonlinear  stage  of  evolution  first  and, 
hence,  bypass  the  evolution  into  Floquet  modes. 

6.  Conclusions  and  discussion.  In  this  paper,  we  have  investigated  the  effect  of  long- 
wavelength  Klebanoff  fluctuations  on  the  instability  of  Blasius  boundary  layer.  By  using  an 
asymptotic  approach  based  on  the  high-Reynolds-number  assumption,  we  derived  a  consistent, 
albeit  simplified,  mathematical  model,  which  appears  to  capture  certain  key  elements  of  this  prob¬ 
lem.  Specifically,  our  analysis  indicates  that  relatively  weak  Klebanoff  fluctuations,  which  do  not 
alter  the  velocity  profile  by  0(1),  may  change  the  near- wall  curvature  of  the  underlying  Blasius  flow 
by  0(1).  This,  in  turn,  has  the  effect  of  modifying  and  even  fundamentally  altering  the  instability 
characteristics  of  the  boundary-layer  flow.  Specifically,  the  perturbed  flow  can  support  primarily 
two-dimensional  instability  modes  that  may  be  localized  in  the  spanwise  direction.  The  spanwise 
distribution  of  these  modes  is  controlled  by  the  distortion  via  a  Schrodinger  equation  governing 
the  mode  shape,  in  which  the  local  skin  friction  of  the  Klebanoff  distortion  acts  as  a  potential. 
The  growth  rates  of  these  modes  are  determined  by  the  wall  torsion  of  the  perturbed  flow. 

When  the  distortion  is  just  strong  enough  to  produce  an  excess  growth  comparable  with  viscous 
growth,  the  abovementioned  instability  modes  may  be  viewed  as  modified  T-S  waves.  However,  as 
the  strength  of  the  distortion  exceeds  a  threshold  range  (in  an  asymptotic  sense),  the  instability 
becomes  essentially  inviscid,  and  the  characteristic  frequencies  and  growth  rates  are  now  much 
higher  than  those  of  the  T-S  waves  in  an  unmodified  Blasius  flow.  A  localized  distortion  (such  as 
what  may  be  induced  when  the  spanwise  correlation  length  of  the  free-stream  disturbance  is  small) 
may  induce  both  sinuous  and  varicose  modes  of  instability.  However,  the  sinuous  modes  are  found 
to  be  more  unstable,  in  general.  Because  the  Klebanoff  distortion  modulates  the  base  flow  in  both  t 
and  x,  its  effect  on  the  instability  is  intermittent  in  time  and  localized  in  space,  i.e.  it  is  manifested 
only  during  a  certain  phase  of  the  modulation  and  in  a  limited  window  along  the  streamwise 
direction.  For  a  localised  Klebanoff  distortion,  for  instance,  the  dominant  sinuous  modes  appear 
only  during  that  phase  in  which  the  distortion  is  characterised  by  a  dominant  low-speed  streak. 
A  spanwise  periodic  distortion  is  found  to  support  analogous,  but  spatially  quasi-periodic  modes 
through  parametric  resonance.  Even  though  we  assumed  the  Klebanoff  distortion  to  have  a  single 
frequency  for  the  most  part,  the  theory  itself  is  quite  general,  being  equally  relevant  to  a  wider 
spectrum  of  vortical  free-stream  disturbances. 

An  interesting  feature  of  the  instability  modes  identified  in  this  paper  is  that,  despite  the  low- 
frequency  nature  of  that  Klebanoff  distortion,  the  unsteadiness  of  the  distortion  plays  a  crucial 
role  in  this  model.  (Leib  et  al.  [24]  had  shown  earlier  that  it  is  also  important  to  account  for 
this  unsteadiness  for  correctly  predicting  the  Klebanoff  distortion  itself.)  Specifically,  the  above 
instability  modes  would  not  have  been  present  for  a  small-amplitude  Klebanoff  distortion  if  it  was 
treated  as  being  steady.  Moreover,  the  intermittent  nature  and  the  convection  of  unstable  patches 
or  spots  are  both  attributed  to  the  unsteadiness  of  the  distortion. 
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The  nonlinear  development  of  a  localised  sinuous  instability  mode  was  also  studied.  In  stage 
I  of  this  evolution,  the  mode  amplitude  is  governed  by  an  integro-or dinary-differential  equation, 
the  solution  of  which  develops  a  finite-distance  singularity  if  the  viscous  effect  is  not  too  large, 
but  decays  if  the  viscosity  parameter  exceeds  a  critical  value.  In  the  vicinity  of  the  singularity, 
the  disturbance  enters  the  second  regime,  in  which  its  evolution  takes  place  over  a  much  shorter 
scale  and  is  governed  by  integro-partial-differential  equation.  While  in  the  first  stage  the  spanwise 
modal  shape  remains  unaffected  by  nonlinearity,  in  the  second  regime,  it  undergoes  nonlinear 
deformation  due  to  the  simultaneous  modulation  by  the  distortion  and  nonlinear  effects.  The 
solution  appears  to  develop  yet  another  singularity  during  nonlinear  stage  II,  which  completes  the 
metamorphosis  of  the  initially  linear  instability  into  fully  nonlinear,  inviscid  disturbances  that  are 
generic  to  the  strongly  nonlinear  phase  of  disturbance  evolution.  It  appears  reasonable  to  suggest 
that  the  temporarily  intermittent  and  spatially  localised  instability  and  its  subsequent  nonlinear 
development  lead  to  patches  of  streak  oscillations  and  even  turbulent  spots.  ^From  our  calculation, 
it  was  estimated  that  the  convection  velocity  of  such  patches  was  about  0.821/oc. 

The  qualitative  predictions  of  our  theory  are  consistent  with  laboratory  observations.  Our 
results  indicate  that  the  streaks  can  become  unstable  even  without  appreciable  changes  in  the 
Blasius  profile.  This  is  precisely  what  Matsubara  et  al.  [25]  concluded  on  the  basis  of  their  experi¬ 
mental  studies.  (See,  also,  Matsubara  &  Alfredsson  [26].)  The  predicted  patches  of  oscillations  and 
turbulent  spots  have  been  observed  in  numerous  experiments.  The  elevated  growth  of  instability 
wavepackets  in  the  presence  of  Klebanoff  distortion  has  also  been  observed  in  the  experiments  by 
Kendall  [20]. 

The  seemingly  puzzling  and  conflicting  experimental  observations  can  be  reconciled  to  some  de¬ 
gree  when  reinterpreted  in  the  light  of  our  theoretical  results.  As  mentioned  in  §1,  the  wavepackets 
develop  out  of  the  background  disturbance  and  amplify  downstream.  Kendall  [19]  associated  these 
with  T-S  waves.  We  believe  that  they  are  likely  to  be  packets  of  the  local  T-S  waves  identified  in 
this  paper  rather  than  the  usual  T-S  waves  in  an  unmodified  Blasius  flow.  These  local  T-S  waves 
exhibit  virtually  all  of  the  unusual  characteristics  as  observed  by  Kendall:  the  onset  threshold, 
the  excess  growth,  and  the  range  of  higher  frequencies.  Since  their  growth  rates  depend  on  the 
magnitude  of  the  Klebanoff  fluctuation,  it  is  to  be  expected  that  their  amplitude  at  a  particular 
observation  point  should  have  a  nonlinear  relation  with  the  magnitude  of  the  Klebanoff  motion. 
Since  the  spanwise  extent  of  these  local  T-S  modes  is  determined  by  the  Klebanoff  distortion,  it  is 
not  surprising  that  their  lateral  spreading  rate  is  much  smaller  than  that  of  a  usual  T-S  wavepacket. 
What  remains  unknown  is  the  receptivity  mechanism  for  the  packets  of  these  local  T-S  waves. 

The  extreme  sensitivity  of  the  boundary-layer  response  to  harmonic  point  excitation  (Watmuff 
[34],  [35])  could  also  be  explained  in  the  light  of  the  present  work.  In  the  presence  of  Klebanoff 
fluctuation,  a  point  excitation  definitely  generates  local  T-S  modes  as  well.  Therefore,  the  general 
response  cannot  be  represented  as  a  summation  of  the  conventional  T-S  waves  only,  as  was  assumed 
in  the  calculation.  This  may  be  the  reason  why  a  meaningful  comparison  was  denied  unless  the 
Klebanoff  fluctuation  is  substantially  reduced. 

The  effect  of  Klebanoff  distortions  on  a  planar  T-S  wave  appears  to  be  rather  controversial 
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(cf.  Boiko  et  al.  [4],  Watmuff  [36]).  In  this  work,  we  did  not  explore  this  phenomenon  in  detail. 
However,  one  might  expect  that  the  continuous  spectrum  of  (3.28)  might  be  involved  in  explaining 
the  observed  deformation  (i.e.  scattering)  of  the  wave  front  by  the  Klebanoff  fluctuation.  This  could 
be  a  topic  for  future  investigation.  It  might  also  be  useful  to  examine  the  connection  between  the 
instability  modes  examined  herein  and  pure  T-S  waves,  i.e.  whether  (and  to  what  degree)  the  above 
modes  actually  emerge  from  existing  T-S  waves  at  appropriate  locations  (and/or  instants)  during 
the  modulation  cycle  and  eventually  degenerate  into  them  at  the  end  of  the  transient  window  of 
instability  (if  they  cannot  reach  nonlinear  amplitudes  by  then).  Such  issues  can  be  exceedingly 
subtle,  as  demonstrated  by  Hall  [17]  for  the  case  of  a  (spatially  homogeneous)  unsteady  Stokes 
layer. 

Jacobs  &  Durbin  [23],  based  on  their  direct  numerical  simulations  of  bypass  transition,  con¬ 
cluded  that  the  streaks  close  to  the  wall  are  stable.  Only  when  the  streaks  lift  up  to  the  edge  of 
the  boundary  layer  to  form  a  ‘backward  jet’  do  they  break  down  into  turbulent  spots.  Our  results 
seem  to  be  in  conflict  with  this  finding.  There  could  be  a  number  of  reasons  for  the  disagreement. 
It  might  be  that  the  present  instability  is  too  weak  and  that  it  is  masked  by  other  more  vigorous 
processes.  Alternatively,  it  is  plausible  that  the  instability  modes  identified  herein  were  not  excited 
in  their  simulations.  Finally,  the  energy  of  the  free-stream  disturbances  in  their  simulation  is  con¬ 
tained  in  a  band  of  rather  high  frequencies  (an  order-of-magnitude  higher  than  typical  frequencies 
of  T-S  waves)  and,  therefore,  the  streaks  are  not  a  linear  response  to  the  low-frequency  compo¬ 
nents.  Rather,  they  are  generated  nonlinearly  by  the  interaction  of  higher-frequency  components. 
The  question  as  to  whether  this  is  the  cause  behind  the  discrepancy  remains  open  at  this  point. 

The  present  work  is,  of  course,  far  from  being  a  complete  or  quantitatively  accurate  descrip¬ 
tion  of  the  problem.  Nonetheless,  the  simplicity  of  the  current  model,  together  with  the  physical 
insights  derived  from  it,  appear  to  justify  the  assumptions  made  herein.  The  theory,  we  believe, 
sheds  useful  light  on  a  very  complex  process  which  has  so  far  eluded  a  first-principles  explana¬ 
tion.  The  primary  shortcoming  of  the  present  theory  corresponds  to  its  neglect  of  the  spanwise 
ellipticity  of  the  Klebanoff  fluctuation.  The  more  general  problem  for  Klebanoff  distortions  with 
an  0(1)  wavelength,  including  the  effects  of  nonlinearity  and  stochasticity,  is  currently  under  in¬ 
vestigation.  Because  no  analytical  techniques  currently  exist  to  deal  with  stochastic  eigenvalue 
problems  with  0(1)  stochastic  parameters,  we  are  pursuing  a  Monte  Carlo  approach  in  an  attempt 
to  quantify  the  expected  growth  rate  of  the  instability  modes  for  a  given  r.m.s.  amplitude  (and 
wavenumber  /frequency  spectrum)  of  the  free-stream  disturbance. 

The  work  of  XW  was  carried  out  during  his  sabbatical  in  2001  at  ICASE,  NASA  Langley  and 
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hospitality  and  financial  support  of  these  institutions.  It  is  a  great  pleasure  to  thank  Prof.  Par  viz 
Moin,  Prof.  Paul  Durbin,  Prof.  W.  C.  Reynolds  and  Prof.  Peter  Bradshaw  for  helpful  discussions 
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